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FOREWORD 


Major  engineering  projects  have  focused  attention  on  problems  created  by 
the  freezing  of  water  in  earth  materials.  These  problems  received 
detailed  attention  during  the  design  of  the  Trans-Alaska  oil  pipeline,  and 
additional  aspects  have  been  encountered  during  the  design  of  the  Alaska 
Natural  Gas  Transportation  System.  They  also  have  received  attention  in 
the  design  and  construction  of  pipelines  in  Siberia  and  Canada.  The 
problems  arise  from  the  tendency  of  ice  to  form  in  lenses  segregated  from 
the  surrounding  soil.  As  these  lenses  grow,  they  are  capable  of 
generating  large  pressures  that  disturb  the  surface  of  the  ground.  The 
term  frost  heaving  is  used  to  describe  this  phenomenon.  Efforts  to  design 
roadways,  airfields,  buildings,  powerlines,  and  other  structures  to 
prevent  damage  from  frost  heaving  date  from  the  early  1900s. 

Because  of  the  seriousness  of  the  problem,  much  time  and  money  have 
been  expended  in  searching  for  methods  of  predicting  and  control  lir^  ice 
segregation  and  frost  heaving.  Unfortunately,  misunderstanding  and 
controversy  have  developed  due  to  inadequate  knowledge  of  the  physical 
processes  involved,  and  due  to  the  widely  differing  backgrounds  of 
scientists  and  engineers  who  have  been  working  on  the  problem. 
Recognizing  this,  the  Committee  on  Permafrost  of  the  Polar  Research  Board 
established  an  Ad  Hoc  Study  Group  on  Ice  Segregation  and  Frost  Heaving  to 
sunmarize  the  current  state  of  knowledge  and  to  identify  areas  of 
controversy  or  uncertainty  that  additional  research  might  resolve.  This 
report  is  the  result.  Its  intent  is  to  provide  a  clearer  understanding  of 
the  general  nature  of  the  field,  of  the  controversial  questions  that 
remain,  and  of  the  research  required  to  resolve  them. 

It  is  the  hope  of  the  Ad  Hoc  Study  Group  and  the  Polar  Research  Board 
that  this  assessment  will  aid  in  the  resolution  of  the  outstanding 
problems  and  thus  assist  the  scientific  and  engineering  community  in 
preventing  and  controlling  the  adverse  consequences  of  frost  effects. 

Charles  R.  Bentley,  Chairman 
Polar  Research  Board 
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EXECUTIVE  soiuaky 


Frost  heaving  caused  by  ice  segregation  occurs  throughout  the  cold  regions 
of  the  world  where  sustained  low  temperatures  are  commonly  encountered. 
Annual  ground  displacements  of  several  inches  with  cyclic  differential 
ground  pressures  of  hundreds  of  pounds  per  square  inch  are  cannon. 
Adverse  effects  are  widespread  and  costly.  Attempts  to  avoid  the  damaging 
effects  of  frost  heaving  often  dominate  design  considerations  and,  thus, 
the  cost  of  highways,  airfields,  and  the  foundations  of  major  structures, 
large  sums  of  money  have  been  spent  during  the  past  two  decades  on  the 
development  of  procedures  and  techniques  designed  to  understand  and 
predict  the  complicated  processes  of  ice  segregation  and  frost  heaving. 
In  addition  to  federal  and  state  agencies,  many  engineering  firms  and 
consortia  of  firms  in  private  industry  are  contributing  to  this  effort. 

During  the  past  five  years  alone,  it  is  estimated  that  in  excess  of 
$50  million  has  been  spent  for  the  design,  construction,  and 
instrumentation  of  experimental  field  sites  at  remote  locations  and  for 
data  collection  and  data  analysis  to  obtain  information  on  ice  segregation 
and  fro6t  heaving  needed  for  the  design  of  the  Alaskan  Natural  Gas 
Transportation  System.  At  least  an  equivalent  amount  of  money  has  been 
expended  on  collecting  geotechnical  data  along  proposed  right  of  ways  and 
in  evaluating  technical  approaches  to  the  prediction  of  frost  heaving 
hazards  for  pipeline  design. 

This  report  and  the  general  recommendations  that  follow  are  Intended 
to  provide  guidance  and  assistance  to  federal  and  state  agencies,  private 
firms,  and  consortia  in  whose  interest  it  is  to  continue  research  in  the 
scientific,  geotechnical,  and  other  engineering  aspects  of  ice  segregation 
and  frost  heaving.  The  Intent  is  to  accelerate  progress  in  attaining  a 
complete  understanding  of  the  phenomena  associated  with  seasonal  and  long¬ 
term  ice  segregation  and  frost  heaving  and  the  application  of  this 
knowledge  in  the  various  planning,  design,  and  maintenance  activities 
associated  with  buildings,  facilities,  airfields,  roads,  utilities,  pipe¬ 
lines,  and  other  structures  wherever  the  need  exists. 

Engineering  hazards  due  to  the  displacements  and  pressures  generated 
by  loe  segregation  and  frost  heaving,  together  with  the  adverse  effects  of 
accumulations  of  segregated  Ice  in  freezing  soil,  are  common  throughout 
most  of  North  America,  Europe,  and  Asia.  Mountainous  areas  that  are 
subject  to  seasonal  freezing  and  thawing  also  are  affected.  Recognition 
of  damage  caused  by  frost  heaving  dates  back  to  Roman  times.  Increasingly 


it  has  been  realized  that,  although  the  engineering  techniques  required  to 
avoid  or  protect  against  frost  heaving  damage  may  be  expensive,  failure  to 
protect  adequately  against  these  effects  often  leads  to  much  higher  costs 
for  subsequent  repair  or  reconstruction. 

Systematic  scientific  investigations  of  frost  heaving  began  in  the 
early  1900s.  Major  advances  in  the  understanding  of  ice  segregation  and 
frost  heaving  were  made  during  the  1950s.  During  the  1970s,  large 
engineering  projects  in  North  America,  Europe,  and  Asia  produced  an 
intensification  of  effort.  As  a  result,  the  principal  processes  of  the 
ice  segregation  that  cause  frost  heaving  have  been  recognized  and  are  now 
being  defined  in  detail.  Freezing  air  temperatures  create  a  thermal 
gradient  that  induces  upward  heat  flow.  Soil  water  freezes  near  the 
ground  surface,  and  segregated  ice  crystals  form  and  grow.  They  then 
coalesce  into  lamelar  lenses  (ice  segregation).  Ice  lens  growth  is  a 
consequence  of  a  continuous  upward  flow  of  water  from  below.  Heat  flow 
and  the  flow  of  soil  moisture  occur  in  the  same  direction  and  are 
interdependent.  During  this  process  a  balance  is  set  up  between  the 
dissipation  of  the  latent  heat  of  freezing  and  the  upward  flow  of  soil 
water.  The  upward  displacement  of  the  soil  surface  (frost  heaving)  is  a 
consequence  of  the  growth  of  a  series  of  Ice  lenses.  When  displacement  is 
restricted,  "heaving  pressures"  develop;  under  certain  conditions,  very 
large  frost  heaving  pressures  are  possible. 

The  process  of  ice  segregation  involves  complex  lnterrelatlon-ahlps 
between  the  ice,  an  unfrozen  liquid  phase,  and  the  bulk  pore  water. 
Complex  latent  heats  of  phase  change  as  well  as  variable  interfacial 
energies  between  phases  are  involved.  Although  there  is  general  agreement 
on  the  existence  of  the  various  phases  and  the  general  concepts  for  which 
governing  equations  and  "rules  of  behavior"  are  being  sought,  we  still  are 
far  from  the  detailed  understanding  needed  to  explain  fully  ice 
segregation  and  to  predict  adequately  rates,  total  and  differential 
displacements,  and  pressures  that  result  from  frost  heaving. 

In  order  to  organize  the  conceptual  information  available  and  to 
utilize  it  for  the  prediction  of  the  various  consequences  of  ice 
segregation  and  frost  heaving,  the  techniques  of  mathematical  simulation 
and  numerical  modeling  have  been  applied.  There  is  widespread  agreement 
that  the  models  of  ice  segregation  and  frost  heaving  must  include 
appropriate  equations  describing  heat  flow  and  the  flow  of  soil  water  and 
that  they  must  deal  with  total  and  differential  heave  displacement  and 
with  frost  heaving  pressures.  Beyond  this  general  agreement,  however, 
there  are  strongly  divergent  views  on  many  fundamental  questions;  for 
example,  the  form  and  basic  structure  of  "an  adequate  model"  is  being 
debated.  There  is  also  disagreement  on  how  to  deal  with  surcharge  and 
overburden  pressures  during  the  ice  segregation  process,  etc.  Finally, 
elaboration  of  existing  one-  and  two-dimensional  models  to  three 
dimensions  is  an  undertaking  of  great  complexity  and  difficulty.  There 
are  many  detailed  questions  to  be  answered  before  significant  progress 
toward  a  complete  mathematical  simulation  of  ice  segregation  and  frost 
heaving  can  be  developed. 


1.  We  recqwnoTH  a  continued  concentrated  research  effort  to 
understand  the  basic  processes  of  ice  segregation  and  frost  heaving. 
Among  the  most  critical  unresolved  questions  and  problems  are: 

•  What  are  the  basic  relationships  that  govern  the  coupled  transport 
processes  (liquid  water,  heat,  water  vapor,  ice)  in  ice  segregation? 

•  How  are  overburden  pressures  related  to  ice  segregation  and  frost 
heaving? 

•  Do  the  same  basic  relationships  between  overburden  pressure  and 
other  parameters  that  characterize  frost  heaving  apply  under  both  static 
and  dynamic  conditions? 

•  What  are  the  important  basic  principles  and  relationships 
governing  long-term,  as  distinct  from  seasonal,  frost  heaving? 

•  What  are  the  maximum  attainable  frost  heaving  pressures  under  long¬ 
term  frost  heaving  conditions  and  how  are  they  best  controlled? 

•  What  are  the  sources  of  the  energy  involved  in  frost  heaving  and 
how  are  they  partitioned  among  the  coupled  processes  involved? 

•  Solute  redistribution  occurs  during  ice  segregation.  More 
information  is  needed  on  the  basic  principles  and  proccesses. 

•  How  is  regelation  involved  in  the  buildup  and  maintenance  of  frost 
heaving  pressures? 

•  What  determines  the  temperatures  and  governs  the  phenomena  of  ice 
nucleation  in  freezing  soils? 

2.  We  reconroend  a  continued  intensive  research  effort  to  develop 
predictive  numerical  models  that,  by  simulating  ice  segregation  and  frost 
heaving,  can  provide  a  means  of  consolidating  all  the  processes  and 
interactions  involved.  This  will  continue  to  aid  in  attaining  a  better 
and  more  complete  understanding  of  frost  heaving  and  for  assessing  its 
impact  in  various  settings  and  circunstances.  Among  the  most  critical 
unresolved  questions  and  problems  pertaining  to  this  aspect  are: 

•  What  are  the  essential  and  most  appropriate  parameters  required? 

•  What  are  the  essential  criteria  for  ice  segregation  and  the 
formation  of  ice  lenses  or  enlarging  ice  masses? 

•  How  should  the  coupled  transport  processes  (heat,  liquid  water, 
water  vapor,  ice)  be  most  effectively  and  appropriately  formulated? 

•  How  can  the  locations  of  enlarging  ice  lenses  be  most  effectively 
specified  and  the  processes  of  "secondary  heaving"  be  better  defined? 

•  Must  the  partition  of  energy  among  the  various  processes  involved 
be  specifically  recognized  and  dealt  with?  If  so,  how? 

•  Are  there  useful  characteristic  frequency  distributions  and  auto- 
correlative  functions  that  can  be  generally  agreed  upon  and  worked  into 
the  various  approaches  to  numerical  simulations  of  frost  heaving? 

•  In  both  seasonal  and  long-tern  simulations,  what  is  the  best 
technique  for  dealing  with  and  scaling  time? 


Progress  In  solving  all  of  these  difficult  and  controversial 
questions  is  crucial  to  the  attainment  of  a  complete,  rational  theory  of 
ice  segregation  and  frost  heaving.  Such  a  theory,  however,  is  needed  for 
the  improved  mathematical  simulations  and  numerical  models  capable  of 
producing  reliable  assessments  of  frost  heaving  hazards  in  a  variety  of 
locations  and  circumstances  where  engineering  design  and  construction  must 
take  place. 

3.  A  common  terminology,  clearly  relating  to  and  bridging  existing 
terminologies  in  the  various  disciplines,  is  needed  to  address  the 
phenomena  of  ice  segregation  and  frost  heaving  and  to  facilitate  a  more 
efficient  and  rational  discussion  of  the  processes,  consequences ,  and 
methods  of  ameliorating  their  adverse  effects.  Therefore,  we  recommend 
that  the  several  national  scientific  and  professional  societies  involved 
give  explicit  attention  to  this  need  and  encourage  efforts  to  achieve  an 
appropriate,  common  terminology. 

4.  Improved  instrumentation  to  measure  the  critical  phenomena 
associated  with  ice  segregation  and  frost  heaving  is  needed.  Most  of  this 
instrumentation  must  be  developed  for  unique  uses  and  circumstances.  We 
recommend  that  this  need  be  explicitly  recognized  by  funding  bodies  and 
that  the  design  and  fabrication  of  needed  instrumentation,  in  addition  to 
the  procurement  of  more  conventional  and  readily  available 
instrumentation,  be  a  prominent  part  of  the  budgets  of  ice  segregation  and 
frost  heaving  research  projects. 

5.  We  recommend  that  field  and  laboratory  studies  of  ice 
segregation  and  frost  heaving  be  better  coordinated  in  order  to  accelerate 
progress  in  the  verification  and  improvement  of  the  various  mathematical 
models  currently  available  and  to  accelerate  the  quantification  of  the 
functions  and  parameters  of  the  various  interrelated  processes.  We 
recomaend  that  this  coordination  be  a  part  of  the  protocol  of  all  major 
research  projects  dealing  with  problems  of  ice  segregation  and  frost 
heaving. 


In  recent  years  professional  societies  have  organized  several 
conferences  to  facilitate  the  exchange  of  knowledge  and  data  on  the 
problems  of  ice  segregation  and  frost  heaving.  However,  much  more  could 
be  done  at  relatively  little  cost  during  the  planning  stages  of  major 
research  projects.  The  needed  coordination  could  be  accomplished  without 
compromising  the  essential  requirements  of  obtaining  and  managing 
propixetary  data.  Indeed,  if  attention  is  given  to  this  recommendation, 
the  value  of  the  proprietary  information  obtained  could  be  vastly 
increased. 


INTRODUCTION 


A  large  fraction  of  the  land  surface  north  of  the  equator  Is  unfrozen  in 
sunnier  but  In  winter  freezes  to  some  depth,  however  slight.  Farther  north 
the  ground  is  normally  frozen  all  the  time  but  thaws  to  some  depth  in 
summer.  When  moist  or  wet  soil  freezes,  it  often  is  accompanied  by 
expansion  due  to  a  process  known  as  frost  heaving,  while  thawing  of 
frozen  soil  often  is  accompanied  by  a  process  known  as  thaw-subsidence. 
Heaving  is  always  associated  with  the  appearance  of  bodies  of  more  or  less 
pure  segregated  ice  within  or  outside  the  soil  while  thaw-subsidence  is 
associated  with  the  melting  of  segregated  ice.  Each  of  these  processes 
has  unwanted  or  damaging  effects  on  engineering  works  and  on  agricultural 
crops  and  enterprises. 

The  appearance  of  segregated  ice,  however,  may  not  be  accompanied  by 
heave  if  there  happens  to  be  a  corresponding  amount  of  Internal 
consolidation  of  the  soil  proper.  Processes  Involved  in  the  formation  of 
segregated  ice  and  the  associated  frost  heave  are  much  more  complicated 
than  generally  is  supposed.  Consequently,  the  effects  are  difficult  to 
predict  with  the  levels  of  confidence  that  an  engineer  requires.  While  a 
number  of  theories  have  been  proposed  to  explain  ice  segregation  and  frost 
heave,  there  is  little  agreement  among  scientists  with  respect  to  which, 
if  any,  of  these  theories  is  physically  consistent  and  sound;  there  is 
uncertainty  as  to  how  to  settle  that  disagreement  and  the  degree  to  which 
all  the  difficult  questions  must  be  settled  in  order  to  develop  reliable 
means  of  predicting  heave  that  will  best  serve  the  most  critical  needs  of 
geotechnical  engineers. 

Ice  segregation  and  frost  heaving  occur  naturally  wherever  the 
climate  is  cold  enough  to  freeze  moist,  fine-grained  soils.  Nearly 
everyone  living  in  the  northern  and  southern  temperate  zones  has 
experienced  the  effects  of  ice  segregation  and  frost  heaving  through  the 
destruction  of  roads  and  highways,  the  displacement  of  foundations,  the 
j aiming  of  doors,  the  misalignment  of  gates,  and  the  cracking  of  masonry. 
Many  people  often  have  simply  and  mistakenly  assumed  that  these  effects 
result  solely  from  the  expansion  of  pore  water  on  freezlig.  When 
confined,  water  can  rupture  pipes,  break  bottles,  and  crack  rocks  as  it 
freezes.  However,  most  of  the  destructive  effects  of  frost  heaving  are 
caused  by  "ice  segregation,"  a  complex  process  that  results  from  the 
peculiar  behavior  of  water  and  other  liquids  as  they  freeze  within  porous 
materials.  In  particular,  water  is  drawn  to  the  freezing  site  from 


elsewhere  by  the  freezing  process  itself.  When  this  water  accumulates  as 
ice,  it  forces  the  soil  apart,  producing  expansion  of  the  external  soil 
boundaries,  as  well  as  internal  consolidation.  The  dynamic  process  of 
ice  segregation  and  the  expansion  resulting  from  freezing  of  the  in  situ 
pore  water,  together,  cause  frost  heaving. 

The  costs  of  frost  damage  in  cold  climates  have  not  been  fully 
determined.  Although  these  costs  are  large,  estimates  have  not  been 
tabulated  to  date,  partly  because  of  the  difficulty  of  recognizing  all  the 
instances  of  frost  effects  and  a  tendency  to  ascribe  many  of  the 
consequences  of  frost  heaving  to  other  causes.  For  example,  the  breakuD 
of  roads  in  spring  is  often  attributed  to  a  combination  of  snow  melt  and 
impeded  drainage.  In  fact,  the  real  cause  is  the  impregnation  of  road 
bases  with  segregated  ice  during  winter.  During  the  spring  thaw,  the  ice 
melts,  weakening  the  base  and  causing  the  pavement  to  fail  when  loaded  by 
traffic.  Even  though  during  the  last  20  years  or  more,  design  procedures 
for  highways,  buildings,  powerlines,  and  other  structures  have  been 
dramatically  improved,  frost  heaving  and  thaw  weakening  are  still  corrmon 
occurrences.  Much  yet  remains  to  be  done  to  improve  procedures  to 
mitigate  frost  action  and  reduce  the  high  costs  of  foundation  construction 
in  cold  climates.  Although  sufficient  understanding  now  exists  to  avoid 
damage  due  to  ice  segregation  and  frost  heaving  in  most  cases,  the  design 
procedures  required  are  often  conservative  and,  therefore,  costly  to 
employ.  Also,  substantial  damage  to  plants  and  crops  occurs  each  year 
wherever  frost  heaving  can  occur.  Further  refinements  in  our  knowledge  of 
ice  segregation  and  frost  heaving  are  needed  to  make  design  criteria  more 
discriminate  and,  thus,  less  expensive  and  to  mitigate  adverse  effects  on 
crops  and  other  agricultural  enterprises. 

Most  of  the  cold  regions  of  the  world  where  the  effects  of  ground 
freezing  are  critical  to  design  considerations  lie  in  the  northern 
hemisphere  continents  of  North  America,  Europe,  U.S.S.R.,  China,  and 
Asia.  The  geographical  areas  sub.ject  to  annual  frost  heaving  and  frost 
effects  are  shown  in  Figure  1.  In  these  areas,  the  design  of  structures 
requires  careful  attention  to  avoid  damage  due  to  frost  action.  The  needs 
of  highway  engineers  have  been  served  by  rules  of  thumb  (based  on  grain- 
size  distributions)  suggested  50  years  ago  by  Casagrande.  The  rules  are 
imperfect  and  surely  result  in  excessively  costly  construction  much  of  the 
time,  but  since  there  can  be  no  assurance  that  adequate  precautions  have 
been  taken  the  rest  of  the  time,  they  are  followed  wherever  possible. 
Casagrande' s  rules  were  devised  for  use  in  projects  that  are  close 
analogues  of  highway  construction.  Special  problems  of  frost  heaving 
induced  by  buried  chilled  pipelines,  for  example,  cannot  be  solved  using 
highway  design  criteria,  yet  such  projects  are  being  planned  for  Alaska 
and  Canada  with  estimated  construction  costs  in  the  range  of  tens  of 
billions  of  dollars.  Clearly  a  much  better  understanding  is  required  for 
large  projects  of  this  type. 

In  the  1950s,  defense  construction  led  to  improvements  in 
geotechnical  engineering  techniques  developed  during  World  War  II  to  cope 
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with  the  Adverse  effects  of  freezing  and  thawing  of  the  ground,  taring 
this  period,  there  was  a  great  demand  for  improved  standards  for  road  and 
runway  construction  on  soils  subject  to  freezing  and  susceptible  to  frost 
heaving,  these  developments  and  the  requirements  for  building  foundations 
on  permafrost  stimulated  research  into  the  rates  and  extent  of  freezing 
and  thawing,  the  dependence  of  these  rates  on  climate  and  ground 
conditions,  and  the  fundamental  principles  governing  the  freezing  and 
thawing  of  porous  media.  The  understanding  gained  during  this  period, 
although  useful,  was  largely  qualitative;  the  engineering  solutions 
devised  as  a  result  were  relatively  simple,  but  they  have  been  effective 
in  many  diverse  circumstances. 

In  the  1970s  a  change  of  emphasis  occurred,  taring  this  period, 
natural  resource  development  in  regions  of  permafrost  led  to  renewed 
attention  to  frost  heaving  phenomena.  The  largest  of  the  projects  of  this 
period,  the  Trans-Alaska  oil  pipeline,  not  only  required  detailed 
techniques  to  avoid  disturbing  the  permafrost,  but  careful  attention  to 
problems  associated  with  frost  heaving  in  the  soil  region  experiencing 
annual  summer  thaw  and  winter  refreezing  and  throughout  extensive  areas  of 
discontinuous  permafrost.  This  attention  continues  in  connection  with  the 
design  of  cold,  burled  natural  gas  pipelines,  burled  tanks  for  the  storage 
of  liquefied  natural  gas,  and  the  artificial  freezing  of  soil  for  purposes 
of  stabilization.  Interest  is  not  confined  to  North  America.  In  northern 
areas  of  the  Soviet  Union,  a  number  of  large-diameter  pipelines  have  been 
constructed  in  the  last  six  years,  and  multiple,  large-diameter  pipes 
running  in  close  proximity  pass  from  the  Urengoy  fields  toward  the  south. 
The  multiple  parallel  lines  provide  a  redundancy  that  is  partly  related  to 
anticipated  difficulties  associated  with  the  possible  thawing  of 
permafrost  and  to  freeze-thaw  phenomena. 

Pipelines  are  particularly  susceptible  to  damage  in  freezing  or 
thawing  ground  because  of  the  extensive  area  and  variety  of  soil  materials 
and  hydrologic  conditions  along  these  routes.  Successful  construction  and 
reliable  operation  require  appropriate  geotechnical  designs  that  deal 
adequately  with  freeze-thaw  phenomena.  In  many  projects,  the  final  costs 
of  pipelines  are  increased  significantly  because  of  the  need  to  deal  with 
the  problems  associated  with  the  freezing  and  thawing  of  the  active  layer 
in  permafrost  and  of  seasonally  frozen  ground  in  areas  where  permafrost  is 
absent. 

In  the  case  of  natural  gas  pipelines  in  permafrost  regions,  plans  may 
Involve  chilling  the  gas  to  below  freezing  temperatures  to  increase 
capacity  and  to  prevent  thawing  of  permafrost.  In  regions  of 
discontinuous  permafrost,  however,  chilling  results  in  potential  danger  of 
frost  heaving  wherever  the  pipe  passes  through  unfrozen  ground.  This 
problem  is  further  compounded  by  the  very  long-term  freezing  period, 
usually  many  years.  In  the  case  of  the  proposed  Alaska  Natural  Gas 
Transportation  System,  as  much  as  20  percent  of  the  alignment  in  Alaska  is 
subject  to  this  hazard. 


Problems  associated  with  frost  heaving  were  not  satisfactorily 
resolved  In  the  design  of  a  proposed  Mackenzie  Valley,  large-diameter  gas 
pipeline  before  this  project  was  cancelled.  During  public  review,  it 
became  apparent  that  the  possible  build-up  of  frost  heave  pressures  due  to 
Ice  segregation  was  being  underestimated.  The  uplifting  forces  that  could 
be  developed  from  the  freezing  of  soil  around  the  pipe  were  initially 
considered  to  be  snail  enough  to  be  counteracted  by  the  weight  of 
overlying  materials.  Additional  Investigations,  however,  demonstrated  the 
possible  development  of  very  large  pressures.  As  a  consequence,  major 
design  changes  including  more  than  surcharge  effects  appeared  to  be 
necessary.  This  matter  still  is  controversial. 

The  design  of  the  Alaska  Natural  Gas  Transportation  System  (known  as 
the  Alaska  Highway  Gas  Pipeline  in  Canada)  has  been  under  development  in 
the  United  States  and  Canada  for  some  time.  During  the  last  five  years,  a 
program  involving  full-scale  test  sites  in  Alaska  has  been  in  progress  to 
develop  appropriate  designs  and  design  procedures  for  the  regions  of 
continuous  and  discontinuous  permafrost.  Expenditures  on  these  test  sites 
probably  have  exceeded  $50  million.  At  the  same  time,  additional  analyses 
involving  mechanical  properties  of  frozen  ground,  as  well  as  heat  flow  and 
moisture  migration,  have  all  been  conducted  at  considerable  expense. 

Another  major  activity  requiring  improved  knowledge  of  the  ice 
segregation  process  is  the  practice  of  artificial  ground  freezing.  This 
is  an  industry  that  has  been  expanding  in  recent  years  as  economic  factors 
have  made  it  competitive  with  other  soil  stabilization  techniques. 
Artificial  ground  freezing  is  used  to  stabilize  soft  soils  and  prevent 
seepage  during  excavation  for  structures,  tunnels,  and  shafts. 
Calculation  of  the  thickness  of  freezewalls  requires  accurate  methods  of 
estimating  moisture  and  heat  flow  during  freezing.  However,  precautions 
must  be  taken  to  prevent  frost  heave  and  thaw  settlement  which  could  cause 
extensive  damage  to  adjacent  structures. 

As  the  behavior  of  natural  soils  during  freezing  is  still  largely 
unknown,  the  need  for  continuing  fundamental  investigations  remains.  A 
basic  requirement  is  the  capacity  to  make  quantitative  predictions  of  the 
magnitude  of  ice  segregation  (Including  the  extent  to  which  this  continues 
in  already-frozen  ground)  and  of  the  pressures  that  can  be  generated  at 
specific  locations  and  in  specific  circumstances.  These  predictions 
require  an  ability  to  predict  accurately  the  flow  of  both  heat  and 
moisture  as  the  ground  freezes  and  thaws  in  a  wide  variety  of  gecmorphic 
and  climatic  conditions. 


CURRENT  KNOWLEDGE 


Early  Investigations 


Frost  heaving  was  first  subjected  to  detailed  study  by  Taber  (1916,  1929, 
1930).  He  observed  that  when  soil  materials  are  frozen  from  the  top  down, 
under  conditions  allowing  water  to  move  upward  within  the  pore  space  from 
below,  the  soil  surface  moved  upward  and  segregated  ice  layers  appeared. 
As  freezing  progressed,  he  could  observe  no  marked  change  in  the  rate  of 
upward  movement  as  one  ice  lens  ceased  to  grow  and  another  one  appeared 
beneath  it.  He  also  noted  similar  behavior  could  be  observed  in  soil 
materials  wetted  with  liquids  such  as  benzine  or  nitrobenzine.  This 
finding  is  important  because  these  liquids  differ  from  water  in  that  they 
contract  in  changing  from  the  liquid  to  the  solid  phase.  This  observation 
establishes  that  ice  segregation  and  the  resultant  frost  heaving  are  due 
principally  to  the  inhibition  and  local  accumulation  of  the  pore  fluid  and 
not  simply  to  the  expansion  of  the  pore  fluid  during  solidification. 

The  principal  processes  of  ice  segregation  as  they  now  are  understood 
are  illustrated  in  Figure  2.  Air  temperatures  below  freezing  create 
freezing  at  the  soil  boundary  and  induce  a  thermal  gradient  in  the  soil 
that  causes  water  flow  toward  the  soil  boundary.  Freezing  is  initiated 
and  the  ice  crystals  grow  to  form  lamelar  masses.  Ice  lens  growth  is 
achieved  by  the  upward  movement  of  water  toward  the  freezing  surface.  A 
balance  is  struck  between  the  rate  of  heat  release  due  to  ice  lens  growth 
and  the  dissipation  by  other  processes  of  the  latent  heat  of  freezing. 
The  displacement  of  soil  toward  the  soil  surface  is  a  consequence  of  ice 
lens  enlargement  and,  when  restricted,  leads  to  the  development  of  heaving 
pressures.  Soil  displacement  occurs  after  soil  consolidation  is  achieved 
and  continues  as  long  as  a  favorable  balance  exists  between  the  rate  of 
heat  removal  and  the  upward  movement  of  soil  water.  Confining  pressures 
and  soil  characteristics  also  are  important.  When  the  heat  and  moisture 
fluxes  became  sufficiently  unequal,  new  ice  lenses  form  at  sites  below, 
where  a  favorable  balance  between  these  two  interacting  fluxes  is 
possible. 

The  ground  surface  depicted  in  Figure  2  shows  a  vertical  cut  to 
emphasize  that  (a)  flow  lines  for  heat  conduction  usually  exit 
perpendicular  to  the  soil  surface,  (b)  flow  lines  for  soil  water  movement 
in  general  are  parallel  to  the  flow  lines  for  heat  conduction,  and  (c)  the 
resulting  planar  ice  lenses  tend  to  be  perpendicular  to  the  soil  water 
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Figure  2.  A  schematic  diagram  of  ice  segregation  and  frost  heaving 
(Anderson  et  al.  1984). 
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flux  and  thermal  flux  flow  lines.  Thus,  it  is  possible  for  the  process  of 
ice  segregation  to  produce  both  vertical  and  horizontal  earth 
displacements.  The  general  rule  is:  displacements  due  to  ice  segregation 
tend  to  be  perpendicular  to  the  ground  surface.  In  general,  the  nearer 
the  water  table  is  to  the  surface,  the  more  readily  water  is  transported 
through  the  partially  saturated,  capillary  zone  to  the  freezing  zone. 

In  addition  to  obtaining  an  early  general  comprehension  of  these 
basic  processes  and  relationships,  Taber  also  recognized  that  the  ice 
lenses,  for  the  most  part,  are  separated  from  the  solid  particles  making 
up  the  soil  structure  during  their  enlargement.  This  is  the  only  way  he 
could  explain  the  clear,  segregated  ice  lenses  that  he  observed  to  form. 
He  reasoned  that  as  heat  was  conducted  toward  the  surface  there  must  be  an 
unfrozen  film  of  water  supporting  the  ice  lenses  as  they  formed, 
separating  them  from  the  soil  fabric.  This  film,  he  visualized, 
maintained  itself  by  inbibing  additional  water  from  the  moist,  unsaturated 
soil  below,  setting  up  a  hydraulic  gradient  that  induced  the  continual 
flow  of  water  required  to  perpetuate  the  process.  Beskow  (1935)  enlarged 
upon  this  basic  concept  and  proposed  that  the  "capillary"  characteristics 
of  soils,  determined  primarily  by  grain-size  distribution,  govern  the 
extent  to  which  soils  may  heave  when  subjected  to  freezing  conditions. 

Some  peculiarities  of  the  process  of  frost  heaving  can  be  illustrated 
by  a  description  of  the  behavior  of  "coarse,"  "medium,"  and  "fine-grained” 
soils  in  either  "open"  or  "closed"  systems  that  may  be  either  "saturated" 
(i.e. ,  virtually  air-free)  or  "unsaturated,"  and  for  cases  in  which  the 
load  supported  by  soil  in  the  zone  of  ice  formation  is  "small"  or  "large." 
Coarse  sand  or  gravel,  a  silt,  and  a  swelling  clay  are  examples  of  what  is 
meant  by  coarse,  medium,  and  fine,  respectively.  Small  loads  may  approach 
zero;  large  loads  may  be  of  the  order  of  1  MRa.  A  closed  system  is  one 
that  is  isolated  from  an  external  source  of  water,  something  that  is  not 
true  of  an  open  system. 

If  loads  are  large  enough,  and  if  the  systems  are  saturated  and 
closed,  all  soils  behave  superficially  in  ways  that,  from  the  engineer's 
viewpoint,  are  practically  the  same.  At  any  stage  of  freezing,  the  heave 
can  be  equated  to  the  volume  change  of  that  part  of  the  initial  water 
present  that  has  changed  to  ice.  In  this  case,  there  may  be  marked 
differences  in  the  appearances  of  the  coarse,  medium,  and  fine  soils,  with 
a  lot  more  ice  segregation  in  the  very  fine  (highly  compressible)  soil 
than  in  the  two  virtually  incompressible  soils  which,  accordingly,  contain 
much  more  "pore  ice." 

If  the  systems  are  open  but  otherwise  as  above,  then  the  coarse  soil 
will  not  heave  at  all;  neither  will  the  medium  soil,  but  the  fine  one  may 
behave  essentially  as  it  did  when  closed,  assuming  that  time  allowed  is 
relatively  small.  If  the  loads  are  very  small,  and  the  systems  are  again 
closed  and  saturated,  the  coarse  soil  will  heave  almost  exactly  as  if  the 
load  was  very  large,  and  so  will  the  very  fine  soil.  The  medium  soil  will 
heave  quite  a  lot  more  than  either  of  the  other  two,  however.  On  the 
other  hand,  if  the  loads  are  small  and  the  systems  are  open  and  saturated 


or  near  saturation,  the  coarse  soil  still  will  not  heave  at  all;  the  clay 
will  behave  much  as  it  always  does,  but  the  raediun  soil  can  produce  truly 
spectacular  heaving! 


Energies  at  the  Interface  and  Unfrozen  Water 


The  earliest  observations  and  insights  on  ice  segregation  have  been 
confirmed,  expanded,  and  elaborated  upon  by  others.  For  example,  the 
"frost  susceptibility"  of  soils  was  extensively  investigated  by  Casagrande 
and  later  by  Kaplar  (1968);  the  role  of  surface  tensions,  particularly  the 
surface  tension  of  the  ice-water  interface,  was  examined  by  Gold  (1957), 
Miller  et  al.  (1960),  Everett  (1961),  Koopman's  and  Miller  (1966),  Miller 
(1973,  1977,  1978),  and  Fenner  (1977),  among  others.  Hie  nature, 
properties,  and  behavior  of  the  transition  zone  at  the  interface 
separating  the  ice  lenses  from  the  soil  mass  also  have  been  the  subjects 
of  many  subsequent  investigations.  For  example,  Bouyoucos  and  McCool 
(1916)  and  Bouyoucos  (1917)  established  that  frozen  soils  retained  a 
fraction  of  their  total  water  content  in  the  unfrozen  state  and  indirectly 
correlated  the  unfrozen  water  with  the  interaction  occurring  at  the 
interface  between  the  mineral  matrix  and  the  soil  water.  Similar 
investigations  were  conducted  by  Beskcw  (1935)  and  Tsytovich  (1945). 

The  early  calorimetric  method  of  determining  the  quantity  of  unfrozt'i 
water  in  frozen  soil  was  refined  by  Martynov  (1956),  Nersesova  aid 
Tsytovich  (1966),  and  Williams  (1963,  1964ab,  1967).  Additional  methods 
of  verifying  the  existence  of  the  unfrozen  water  in  capillary  spaces  and 
Interfacial  films  and  of  defining  the  properties  of  the  unfrozen  liquid 
water  in  these  domains  were  devised  by  Anderson  and  Hoekstra  (1965ab), 
Anderson  (1966,  1967a),  and  Anderson  and  Tice  (1972).  Less  direct  methods 
that  corroborate  the  existence  of  the  unfrozen  water  have  been  contributed 
by  Hoekstra  (1965),  Hoekstra  and  Chamberlain  (1964),  Nakano  et  al.  (1972), 
and  Riitterson  and  Sknith  (1980). 

As  a  result  of  these  investigations,  the  liquid-like  characteristics 
of  the  water  in  frozen  soil  materials  have  been  described  and  well 
established.  Electrical  conductance  measurements  have  demonstrated  the 
mobility  of  ionic  species  in  frozen  soil  (Hoekstra  1965,  1909a;  Murrmann 
et  al.  1970).  The  fact  that  water  in  frozen  soil  can  be  mobilized  by 
electrical  or  thermal  gradients  demonstrates  the  continuity  and  the  high 
molecular  mobility  of  the  water  in  the  Interface  between  the  ice  and  the 
soil  particles.  That  soil  particles  can  migrate  through  ice  masses  proves 
that  Individual  mineral  particles  are  free  of  any  but  the  most  feeble  and 
transitory  connections  to  the  ice  phase  (Hoekstra  1969a).  Nuclear 
magnetic  resonance  (NMR)  data  indicate  that  the  transition  zone  retains 
its  liquid-like  characteristics  to  very  low  temperatures  (-40°C  or 
lower).  The  MR  data  permit  an  estimate  of  the  viscosity  of  the  unfrozen 


water  at  the  Interface.  Although  somewhat  unrealistic  because  of  the 
small  physical  dimensions  involved,  it  has  been  suggested  that  the 
viscosity  of  the  unfrozen  water  is  about  600  times  greater  than  that  of 
ordinary  water,  comparable  to  the  viscosity  of  glycerol  (Fripiat  et  al. 
1965;  Anderson  and  Morgenstem  1973).  Measured  diffusion  coefficients  in 
frozen  soils  are  consistent  with  this  view  (Murrraann  et  al.  1970).  This 
conclusion  is  not  fully  established,  however. 

Typical  relationships  between  the  quantity  of  unfrozen  water  in 
frozen  soil  and  temperature  are  shown  for  three  representative  soil 
materials  in  Figure  3  (Tice  et  al.  1978).  The  figure  shows  that  the 
effect  of  varying  total  water  and  ice  contents  is  minor.  The  quantity  of 
unfrozen  water  present  depends  principally  on  the  temperature.  When  the 
data  of  Figure  3  are  normalized  by  dividing  the  unfrozen  water  contents  at 
any  temperature  by  the  specific  surface  area  of  the  soil  materials,  the 
data  tend  to  coalesce.  Calculated  thicknesses  of  the  unfrozen  water  films 
resulting  from  this  treatment  range  from  approximately  two  diameters  of 
the  water  molecule  at  temperatures  of  -10°C  and  lower  to  progressively 
larger  dimensions  at  temperatures  approaching  the  melting  point.  Other 
factors,  such  as  salinity  and  pressure  also  act  to  influence  the  thickness 
of  these  unfrozen  water  films  and  the  quantity  of  unfrozen  water  present 
at  any  given  temperature.  An  increase  in  salinity  or  pressure  increases 
the  quantity  of  unfrozen  water  at  a  given  temperature  (Banin  and  Anderson 
1974;  Anderson  1967b). 

It  is  possible  to  visualize  a  variety  of  physical  arrangements  of  the 
interconnected,  unfrozen  water  films  throughout  a  frozen  soil  mass.  The 
experimental  investigation  and  analysis  of  Koopmans  and  Hiller  (1966) 
relate  both  temperature  and  pressure  to  unfrozen  water  content  in 
saturated  soils.  This  work  suggests  that  under  many  common  conditions 
near  the  site  of  active  segregation  a  substantial  amount  of  pore  water 
remains  unfrozen  in  the  smaller  pores,  in  addition  to  that  in  adsorbed 
films.  Simulations  based  on  the  concepts  of  Hiller  (1978)  support  this 
view  (O'Neill  and  Hiller,  1962,  1964).  This  means  that  unfrozen  pore 
space,  as  well  as  the  unfrozen  films,  may  contribute  significantly  to  the 
pathways  that  liquid  water  may  take  in  feeding  ice  lens  growth. 


Hydrological  Aspects 


The  role  that  ice  lenses  play  in  restricting  infiltration  during  the 
winter  and  spring  melt  is  important  in  governing  surface  water  runoff  and 
groundwater  recharge,  later  from  the  snowpack  can  either  return  to  the 
atmosphere,  infiltrate  into  subsurface  systems,  or  leave  as  surface 
runoff.  The  infiltration  rate  into  seasonally  frozen  soils  is  controlled 
by  the  fraction  of  the  pore  space  that  is  occupied  by  ice  and  by  the 
barrier  to  flow  that  is  created  by  ice  lenses.  Seasonally  frozen  soils 
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Pimire  3  A  plot  of  unfrozen  water  content  versus  temperature  for  several 
representative  imterials  (Andersland  and  Anderson  1978). 
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with  low  moisture  contents  have  relatively  high  infiltration  rates  because 
there  is  little  ice  blocking  the  soil  pores.  This  facilitates  groundwater 
recharge.  Soils  containing  large  quantities  of  ice  are  less  permeable;  in 
this  case  surface  water  runoff  is  enhanced. 

It  is  a  generally  accepted  view  that  ice  lens  development  occurs 
somewhat  behind  the  freezing  front.  As  can  be  visualized  from  Figure  2, 
when  unfrozen  water  is  drawn  upward  to  the  freezing  front  in  response  to  a 
hydraulic  gradient,  there  may  be  a  distinct  zone  between  the  freezing 
front  and  the  region  where  ice  lenses  develop;  this  is  referred  to  as  the 
"frozen  fringe."  The  hydraulic  conductivity  of  the  frozen  fringe  is  one  of 
the  principal  governing  factors  that  determines  where  ice  lens  growth 
actually  occurs.  In  addition  to  a  dependence  on  the  hydraulic 
conductivity  of  the  frozen  fringe,  Konrad  and  Morgenstern  (1981)  have 
related  the  position  where  ice  lens  development  occurs  to  the  temperature 
gradient.  Penner  (1982)  has  shown  that  the  position  of  the  ice  front  is 
also  related  to  the  overburden  pressure;  Increasing  the  overburden 
pressure  lowers  the  temperature  at  the  ice  lens. 

The  hydraulic  gradient  is  composed  principally  of  a  gravitational 
component  and  a  "soil  water  potential"  component.  The  gravitational 
potential  gradient  is  constant,  is  independent  of  the  freezing  process, 
and  is  always  directed  downward.  The  soil  water  potential  gradient,  on 
the  other  hand,  depends  on  the  freezing  process  and  may  be  orders  of 
magnitude  greater  than  the  gravitational  component.  Other  potential 
energy  gradients  that  may  be  involved,  such  as  that  arising  from  osmotic 
forces,  are  important  when  solutes  are  present.  All  things  considered, 
however,  net  movement  of  water  within  the  soil  generally  is  from  warmer  to 
colder  areas.  Thus,  the  direction  of  movement  of  soil  moisture  normally 
coincides  with  the  direction  of  heat  transfer. 

The  process  of  lens  initiation  and  growth  within  the  frozen  soil 
behind  a  frozen  fringe  has  been  referred  to,  or  designated  as,  "secondary 
heave"  (Miller  1972).  The  possibility  of  continual  redistribution  of 
water  among  the  ice  lenses  also  is  an  important  facet  of  frost  heaving 
phenomenology,  later  can  migrate  from  ice  crystals  and  ice  lenses  within 
the  frozen  zone  from  frozen,  but  warmer  parts  of  the  soil,  to  colder 
areas.  This  redistribution  of  soil  water  within  the  frozen  zone  must 
proceed  at  a  much  slower  rate  than  the  upward  movement  of  water  to  the 
freezing  front,  however,  because  of  the  much  lower  hydraulic 
conductivities  of  the  frozen  material.  The  hydraulic  conductivity  of 
saturated  frozen  soils  is  much  lower  than  that  of  the  same  soils 
unfrozen.  As  the  temperature  falls,  the  hydraulic  conductivity  decreases 
because  of  the  freezing  of  capillary  water  and  the  decrease  in  thickness 
of  the  unfrozen  water  films.  The  hydraulic  conductivity  does  not  fall  to 
zero,  however.  Seme  permeability  to  water  remains,  even  at  quite  low 
temperatures. 

The  efficiency  of  ice  segregation  varies  with  soil  texture.  Silty 
soils  generally  account  for  most  engineering  problems;  they  are  the  most 
highly  susceptible  to  ice  segregation.  Under  comparable  conditions,  soils 
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finer  than  silts  generally  contain  less  segregated  ice  and  are  less 
susceptible  to  frost  heaving,  probably  because  their  hydraulic 
conductivities  are  so  much  lover,  both  in  the  saturated  and  the 
unsaturated  condition,  than  those  of  coarser  soils.  The  hydraulic 
conductivity  of  a  given  freezing  soil  ranges  between  the  lov  values 
observed  during  a  freezing  cycle  and  the  high  values  observed  during  a 
thawing  cycle.  Hysteresis  effects  on  all  soil  properties  are  amplified  by 
freezing  and  thawing.  Thus,  the  differences  between  wetting  and  drying 
cycles  are  amplified,  as  are  disparities  between  laboratory  samples. 
Variations  in  the  hydraulic  and  thermal  properties  of  soils,  coupled  with 
variations  in  available  soil  moisture,  result  in  variations  in  the  depth 
of  frost  penetration  and  the  amount  of  frost  heave.  The  differential 
rates  and  amounts  of  frost  heave  from  one  adjacent  location  to  another 
often  are  troubling  aspects  of  frost  heaving.  In  fact,  these  differences 
are  probably  the  principal  cause  of  most  geotechnical  problems  resulting 
from  frost  heaving. 

Heaving  always  results  in  irreversible  changes  in  the  structure  of 
soils.  Mechanical  properties  are  greatly  affected  when  ice  lenses  develop 
and  the  redistribution  of  soil  water  occurs.  Compressible  soils  such  as 
clays  consolidate  beneath  a  growing  ice  lens  if  the  effective  stress 
increase  caused  by  the  pore  water  suction  exceeds  the  past  maximum 
effective  stress.  This  may  not  cause  a  thaw  weakening  problem  during 
subsequent  thawing,  if  melt  water  can  be  readily  drained  away.  In  fact, 
it  is  sometimes  possible  to  take  advantage  of  this  process  to  improve 
strength  and  consolidation  characteristics.  However,  if  the  possibility 
of  draining  during  thaw  is  not  recognized,  estimates  of  overconsolidation 
by  freezing  can  be  highly  inaccurate  (Chantoerlain  and  Blouln  1978).  This 
can  be  a  problem  in  applications  of  artificial  ground-freezing  techniques, 
where  a  common,  erroneous  assumption  is  that  an  estimate  of  thaw 
consolidation  can  be  taken  as  equal  to  measured  frost  heave. 

Ice  segregation  also  causes  changes  in  the  hydraulic  conductivities 
of  soils  in  another  way.  The  structural  changes  caused  by  the  increased 
effective  stress  within  the  "frozen  fringe"  or  in  the  unfrozen  zone 
beneath  the  freezing  front,  are  known  to  cause  large  Increases  in  the 
thawed  hydraulic  conductivities  of  soft  clay  soils,  whereas  the 
disintegration  of  soil  aggregates  that  is  caused  by  ice  segregation  in 
some  silty  soils  results  in  lower  hydraulic  conductivities.  Thus,  ice 
segregation  and  frost  heaving  during  repeated  freeze-thaw  cycles  cause 
changes  in  soil  structure  that  also  affect  hydraulic  conductivity. 

Most  studies  on  frost  heave  have  been  concerned  primarily  with  near¬ 
surface  conditions  and  relatively  short  time  periods.  This  is  appropriate 
in  the  case  of  the  seasonal  frost  heave  of  roads  and  airfields,  and  in  the 
formation  of  needle  ice  in  the  vegetation  mat  and  the  soil  surface.  The 
potential  of  damage  to  a  buried  chilled  gas  pipeline  by  continual  frost 
heave  during  the  service  life  or  to  a  burled  liquefied  gas  tank  is  a  more 
recent  concern  involving  longer  times  and  greater  depths.  This  new 
perspective  has  served  to  broaden  and  extend  considerations  of  frost 
heaving  phenomena. 
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Massive  Ground  Ice 


Excess  ice  produced  by  ice  segregation  is  by  no  means  confined  to  small 
ice  lenses  at  shallow  depths.  There  is  abundant  evidence  from  natural 
exposures  along  river  banks  and  coastal  bluffs,  drilling,  resistivity  and 
gravity  surveys  (Hampton  and  Walcott  1974)  to  show  that  large  tabular 
bodies  of  ice,  commonly  referred  to  as  massive  ice,  are  widespread  in  some 
permafrost  areas  of  North  America  and  the  Soviet  Union.  Most  of  the 
massive  ice  so  far  encountered  in  drill  holes  lies  within  the  5-  to  40-m- 
depth  range  (Mackay  1971,  1973ab),  but  some  massive  ice  has  been  found  at 
depths  exceeding  100  m  (Dubikow  1982).  The  area  of  some  discrete  bodies 
of  massive  ice  may  exceed  0.5  km2. 

The  thickness  of  massive  ice  accumulations  comaonly  ranges  from  less 
than  1  m  to  a  maximum  reported  thickness  of  about  40  m.  Massive  ice 
bodies  show  every  conceivable  gradation  in  purity,  from  ice  with  no 
sediment  at  one  extreme  to  high  ice  content  sediments  at  the  other 
extreme.  In  general,  massive  ice  occurs  near  the  basal  portion  of  fine¬ 
grained  material  that  typically  is  underlain  by  sand  (Mackay  1973ab; 
Dubikow  1982).  The  fine-grained  sediments  in  which  the  massive  ice  is 
found  show  a  wide  range  of  deposltlonal  environments  such  as  basal  till, 
lake  clays,  and  marine  clays.  The  massive  ice  is  usually  layered  or 
stratified  when  viewed  in  natural  exposures  as  if  the  ice  had  grown  in 
situ  and  distended  the  sediments,  in  accord  ism- like  fashion.  Some  of  the 
massive  ice  is  deformed.  Folding  can  often  be  linked  with  differential 
loading. 

There  has  been  considerable  debate  as  to  whether  the  massive  ice  is 
burled  glacier  ice,  segregated  ice  which  has  grown  in  place,  or  ice  frozen 
from  bulk  water  injected  under  pressure  into  frozen  or  unfrozen  sediments 
and  then  frozen  (e.g. ,  Mackay  1971,  1973ab  ;  Karpov  1981;  Kbreisha  et  al. 
1981;  Dubikow  1982;  Anisimova  and  Krltsuk  1983).  Some  massive  ice  is 
doubtless  burled  glacier  ice.  Indeed  some  exposures  of  debris-rich 
regelation  ice  known  to  have  come  from  contemporary  glaciers  closely 
resemble  massive  ground  ice.  However,  much  of  the  massive  ice  probably  is 
not  burled  glacier  ice  because  (a)  some  massive  ice  is  in  unglaclated 
areas,  (b)  some  massive  ice  is  found  in  lacustrine  and  marine  sediments  or 
lies  beneath  topographic  highs  where  burial  would  be  difficult  to  explain, 
(c)  the  salinity  of  most  of  the  massive  ice  is  far  greater  than  that  of 
glacier  ice,  (d)  the  ice  petrofabrlcs  are  similar  to  that  of  thin  lenses 
of  segregated  ice,  and  (e)  vertical  profiles  of  the  ionic  content  of  the 
ice  (water)  from  the  ground  surface  to  the  sediments  beneath  the  massive 
ice  tend  to  show  a  cannon  water  origin  for  the  ice  (water),  all  of  this  is 
incompatible  with  a  glacial  origin. 

Most  massive  ground  ice  probably  is  segregated  ice  frozen  in  an  open 
system  with  the  source  water  being  under  pressure  (e.g.,  Karpov  1981). 
The  common  occurrence  of  massive  ice  in  fine-grained  sediments  underlain 
by  sands  cannot  be  coincidental;  it  suggests  rather  that  the  swnrfa  provide 


the  open  system  for  water  under  pressure  needed  to  sustain  segregated  ice 
growth  for  a  long  period  of  time.  Evidence  for  water  under  pressure  can 
sometimes  be  seen  in  injection  forms  of  pure  ice  such  as  dikes,  sills,  and 
Irregularly  shaped  finger-like  ice  masses.  The  water  quality  analyses  are 
in  agreement  with  what  would  be  expected  from  the  growth  of  segregated  ice 
by  downward  freezing  (Anisimova  and  Kritsuk  1983;  Dubikcw  1982).  The 
mineralization  often  increases  with  depth,  presumably  from  solute 
rejection  during  freezing.  Mineralization  is  also  related  to  soil  type 
and  ice  content.  Oxygen  isotope  analyses  have  been  used  to  obtain  an 
indication  of  the  water  source  of  some  of  the  massive  ice  (Mackay 
1983ab).  Data  available  for  the  Beaufort  Sea  coast  indicate  a  cold 
freshwater  source. 


Mathematical  Simulations  and  Numerical  Modeling 


As  mentioned  earlier,  in  ice  segregation  and  frost  heaving,  transient 
processes  of  heat  and  mass  transfer  are  coupled.  Consequently,  nearly  all 
mathematical  simulations  contain  a  set  of  coupled  partial  differential 
equations  describing  the  flow  of  heat  and  water.  Although  these  equations 
are  inherently  nonlinear,  during  the  course  of  the  various  mathematical 
treatments  employed  they  typically  become  linear.  The  complexity  of  ice 
segregation  and  frost  heaving  requires  numerical  solutions,  typically  with 
finite  differences  applied  in  time,  and  either  finite  differences  or 
finite  elements  applied  in  space.  Hysteresis  complicates  the  problem 
greatly.  In  fact,  no  heave  simulations  are  yet  known  that  are  capable  of 
dealing  with  hysteresis,  although  as  already  pointed  out,  hysteresis  is 
known  to  be  a  conmon  occurrence. 

"Modeling"  is  a  term  that  has  become  fashionable  during  the  last  two 
decades.  "Mathematical  modeling"  refers  to  the  development  of 
mathematical  representations  based  on  physical  concepts,  instead,  of  the 
development  and  construction  of  physical  models.  With  the  advent  and 
present  widespread  availability  of  computers,  mathematical  modeling  has 
become  increasingly  useful. 

The  objective  of  the  modeler  determines  to  a  large  extent  the  type  of 
model  developed.  For  example,  engineering  models  tend  to  be  primarily 
predictive;  models  developed  for  scientific  purposes  tend  to  be  much  more 
elaborate  and  to  reflect  an  effort  toward  a  complete  understanding  of  key 
aspects  of  the  natural  processes  Involved.  Scientists  connonly  take  a 
constructionist  approach  and  attempt  to  model  every  aspect  of  a  process. 
Engineers  often  employ  a  systems  approach  that  can  include  rather  complete 
descriptions  of  certain  physical  and  thermodynamic  aspects  of  the  system 
of  processes,  but  in  the  absence  of  complete  understanding  they  ultimately 
must  resort  to  simpler,  phenomenological  descriptions  of  the  processes. 
In  some  modeling  efforts  only  the  sign  or  order  of  magnitude  change  in  an 
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output  variable  Initiated  by  a  parameterization  change  or  alteration  of  an 
input  variable  or  boundary  condition  is  needed.  This  is  a  kind  of 
"sensitivity  testing"  that  is  frequently  employed  when  the  system 
parameterization  is  weak  due  to  inadequate  Information.  In  other  cases 
detailed  models  are  tailored  for  "site  specific"  design  analysis. 

It  is  possible  for  models  to  produce  simulations  that  overreach 
understanding  and  invite  self-delusion.  Limitations  in  understanding  the 
underlying  processes  and  the  great  variety  of  circumstances  under  which 
frost  heave  occurs  can  lead  to  modeling  excesses  in  addition  to  errors 
from  other  sources.  It  is  also  possible  to  construct  models  that  are  so 
elaborate  that  the  parameters  required  cannot  be  adequately  identified, 
specified,  or  measured.  So  many  Interacting  processes  may  be  Included 
that  cumulative  errors  of  measurement  and  parameter  specification, 
together  with  the  practical  limitations  of  testing  and  application,  lead 
to  grossly  inaccurate  results.  Clear-cut  conclusions  may  be  impossible  to 
extract  from  the  results.  Yet  when  the  physical  complexity  is  avoided  by 
a  series  of  oversimplifications,  similar  difficulties  arise. 

Notwithstanding  these  limitations,  there  is  much  promise  in  the 
modeling  of  ice  segregation  and  frost  heave.  Future  new  models  may  become 
more  complex  in  order  to  achieve  closer  fidelity  to  the  processes  and 
conditions  involved.  In  any  case,  rapid  progress  is  expected  as  the 
limitations  and  pitfalls  are  clearly  recognized  at  each  step  and  if 
applications  are  carefully  chosen. 

All  frost  heave  models  must  handle  the  latent  heat  effects  that 
accompany  phase  changes.  These  effects  are  extraordinarily  large,  and 
they  can  be  troublesome.  Crank  (1981)  provides  a  review  of  some  practical 
methods  for  handling  a  simple  phase  change.  If  one  assumes  that  the  phase 
transition  occurs  as  a  step  change  in  space,  several  methods  recently 
developed  to  attack  Stefan-type  problems  are  useful,  (Bonnerot  and  Jamet 
1979;  Crank  and  Prahle  1973;  lynch  and  O'Neill  1981).  These  methods 
faithfully  model  the  phase  discontinuity  by  ensuring  that  a  numerical  mesh 
boundary  coincides  with  it.  It  is  difficult  to  use  approximate  techniques 
and  finite  mesh  spaclngs  to  deal  with  discontinuities,  and  even  more 
difficult  to  deal  with  moving  discontinuities.  These  techniques  avoid 
approximation  across  the  zone  of  highly  concentrated  activity  and 
discontinuity  and  express  latent  heat  effects  through  boundary  or 
compatibility  conditions. 

Problems  also  may  be  approached  by  considering  the  occurrence  of  the 
phase  transition  as  a  step  in  macroscopic  space  but  amenable  to  a  smoothed 
numerical  representation.  However,  some  systematic  artificial  distention 
of  the  transition  zone  must  be  used  such  that  energy  is  conserved.  This 
approach  is  the  essence  of  a  variety  of  "enthalpy  approaches"  (Voller  and 
Cross  1981).  In  finite  element  techniques  latent  heat  effects  may  be 
spread  over  an  entire  finite  element.  Thus,  the  phase  change  locale  is 
arbitrarily  expanded  to  the  model's  scale  of  discrimination.  This  feature 
is  similar  in  concept  to  some  other  related  techniques  such  as  the  method 
of  "excess  degrees"  of  Dusinberre  (1949)  and  the  approaches  of  Doherty 
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(1970)  and  Hromadka  et  al.  (I960,  1982).  These  approaches  Involve  methods 
of  repeatedly  forcing  computed  freezing-zone  temperatures  back  to  a  phase- 
change  temperature  as  heat  is  withdrawn  until  the  latent  heat  of  that  zone 
is  exhausted.  For  situations  in  which  the  theory  applied  or  the  scale  of 
analysis  requires  or  involves  a  physically  extended  phase-change  zone,  a 
nonstep  freezing-zone  representation  follows  naturally  (e.g.,  Taylor  and 
Luthin  1978;  O’Neill  and  Miller  1982;  Miller  1981). 

The  principle  of  conservation  of  momentum  (balance  of  stress)  also 
must  be  incorporated  in  any  frost  heave  model.  Experimental  evidence 
(e.g.,  Renner  and  Ueda  1979)  shows  that  an  Increase  in  overburden  pressure 
(accumulated  weight  on  the  freezing  zone)  diminishes  the  amount  of  heave 
produced  by  forces  acting  in  the  freezing  zone,  other  things  being  equal. 
Overburden  pressure  is  defined  as  a  macroscopic  net  force  per  unit  area 
that  must  be  balanced  by  the  combination  of  forces  developed  in  each 
constituent  of  the  freezing  or  thawing  soil.  If  a  weight  is  applied  at 
the  top  of  a  soil  column  of  unit  area,  the  overburden  pressure  at  the  top 
equals  this  weight.  At  the  bottom  of  the  column  the  overburden  pressure 
equals  that  weight  plus  the  weight  of  the  column. 

The  pressure  in  any  phase  or  constituent  usually  differs  from  the 
overburden  pressure.  For  example,  if  the  pore  water  pressure  in  a 
saturated,  coarse-grained,  unfrozen  soil  is  zero  gage,  then  any  weight  on 
the  soil  and  the  soil'  column  above  must  be  supported  entirely  by  stresses 
transmitted  through  intergranular  contacts  of  the  soil  mineral.  This  is 
commonly  known  as  the  "effective  stress."  Microscopically,  the  stresses 
caused  by  intergranular  contacts  are  quite  complex.  The  details  may  be 
smoothed,  however,  and  referred  to  simply  as  net  intergranular  force  per 
unit  area  of  soil.  Usually,  the  overburden  pressure  is  not  fully  balanced 
by  the  effective  stress;  pore  contents  participate  as  well.  The 
macroscopic  partitioning  of  pore  stress,  or  so-called  "neutral  stress," 
into  degrees  of  participation  by  the  various  pore  constituents  Involves  a 
microscopically  complex  interplay  of  stresses  within  the  constituents,  of 
energies  at  the  interface,  and  of  microscopic  geometric  configuration. 
Miller  (1978)  has  suggested  the  use  of  a  "stress  partitioning  factor"  to 
weight  the  relative  participations  of  macroscopic  pore  constituent 
stresses.  Recent  work  (Snyder  1980)  has  provided  some  relevant 
experimental  results,  and  this  concept  is  an  integral  part  of  the  model  of 
O'Neill  and  Miller  (1982,  1984). 

Reasoning  by  analogy  with  unfrozen,  swelling  clays,  Groenevelt  and 
Kay  (1977)  have  derived  equations  connecting  overburden  pressure,  ice 
pressure,  liquid  water  pressure,  and  the  volume  ratios  of  ice  and  liquid. 
Sheppard  et  al.  (1978)  refer  to  this  approach  in  the  presentation  of  their 
model  but  do  not  use  it.  Groenevelt  and  Kay  (1980ab)  have  developed 
additional  relationships  between  envelope,  liquid,  ice,  and  solid 
pressures  from  thermodynamic  concepts. 

One  may  question  the  validity  of  construing  stresses  within  the 
system  as  having  simple  macroscopic  pressure-like  manifestations, 
especially  when  the  stresses  refer  to  an  ensemble  of  grains  or  correspond 
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to  forces  within  extremely  thin  adsorbed  layers,  tightly  bound  to  the 
mineral  particles  (Takagl  1980;  Miller  1980ab).  it  has  been  suggested 
that  when  water  In  the  films  pulls  nearby  water  in  while  simultaneously 
pushing  away  the  ice  that  forms,  it  exhibits  a  state  of  stress  which  is 
nonisotropic  and  otherwise  complex  (Miller  1980b;  Vignes  and  Dijkema  1974; 
Vignes-Adler  1977).  Derjaginn  and  Churaev  (1978)  also  invoke  surface 
layer  forces  to  explain  the  expulsion  of  particles  from  adjoining  ice. 
Romkens  and  Miller  (1973)  propose  other  descriptions  of  transfer  processes 
in  the  film  surrounding  mineral  particles  embedded  in  ice.  Osmotic 
pressure  concepts  also  have  been  used  (Miller  1980a),  to  explain  film- 
water/pore-water/  ice-pressure  relationships.  Much  remains  to  be  done 
before  consensus  is  achieved.  Consequently,  defining,  measuring,  and 
utilizing  the  "right"  pressures  remain  formidable  obstacles. 

The  presence  of  unfrozen  water  in  a  heaving  soil  is  taken  into 
account  in  different  ways  by  the  various  frost  heave  models  now  being 
developed.  A3  Indicated  earlier,  a  major  part  of  the  unfrozen  water  flow 
takes  place  within  the  warmest,  least  frozen  portion  of  the  frozen  soil. 
Ice  segregation  usually  occurs  within  the  frozen  zone,  slightly  above  a 
"frozen  fringe."  The  frozen  fringe  is  defined  as  the  zone  in  which 
volumetric  ice  content  increases  from  zero  at  the  freezir^  front  to  nearly 
100  percent  at  the  location  of  the  warmest  lens  (e.g.,  Hoekstra  1969ab; 
Loch  and  Kay  1978).  The  thickness  of  the  frozen  fringe  is  variable  and 
depends  on  several  factors.  Two  are  particularly  Important:  high 

overburden  pressures  and  low  temperature  gradients  both  act  to  increase 
the  thickness  of  the  unfrozen  fringe. 

A  well-known  relationship  that  has  been  accepted  as  applicable  to  the 
ice-liquid  mixture  in  frozen  soil  is  the  generalized  Clapeyron  equation 

Vw  d(Pw  "  ")  *  vi  ^i  "  T  dT  (D 

where  ° 

L  *  latent  heat  of  fusion  of  water  (80  cal/g); 
vi(Vw)  ■  specific  volume  of  ice  (liquid  water)  (m3/kg); 

Pi(Pw)  *  pressure  acting  on  the  ice  (liquid  water)  (N/m2); 

T  •  temperature  (degrees  Kelvin); 

T0  *  freezing  point  of  bulk  water  at  atmospheric  pressure 
(degrees  Kelvin);  and 

»  -  osmotic  pressure  (N/m2). 

Recent  derivations  of  this  equation  have  been  offered  by  Kay  and 
Groenevelt  (1974)  and  Loch  and  Kay  (1978).  This  equation  can  be  used  to 
represent  the  equilibrium  between  pore  ice  and  unfrozen  water  distributed 


between  two  limiting  conditions:  at  one  extreme,  water  in  thin  films  at 
the  Interface,  with  high  overall  pressure  and  high  osnotic  pressure;  at 
the  other  extreme,  unfrozen  pore  water  with  much  lower  overall  pressure 
and  osmotic  pressure  (Miller  198Qab).  When  applied  macroscopically,  Pw 
and  *  in  this  equation  are  regarded  as  pertaining  primarily  to  pore 
water.  The  difference  (Pw  -  ff)  is  regarded  as  unique  for  the  medium 
in  a  given  state. 

A  number  of  investigators  have  noted  the  experimental  evidence 
linking  this  concept  with  capillary  theory.  When  liquid  and  another  water 
phase  coexist  in  a  configuration  with  curved  phase  interfaces  (as  in  soil 
pores),  there  is  a  pressure  difference  between  the  two  phases  due  to 
effects  at  the  interface,  expressible  in  general  as 

Pj  “  *  wjwF(O)  (2) 

where  j  is  either  "i"  for  ice  in  a  saturated  frozen  soil,  or  is  "a"  for 
air  in  an  unsaturated  unfrozen  soil,  and  0  denotes  unfrozen  water 
content.  The  factor  wj*  is  the  surface  tension  (N/m)  associated  with 
the  interface  in  question.  F(0)  is  a  function  that  appropriately  repre¬ 
sents  the  mean  curvature  of  the  microscopic  phase  interfaces. 

The  earliest  applications  of  this  approach  treat  F(0)  as  a  constant. 
Some  of  the  simpler  more  recent  models  still  do  so.  For  example,  one  may 
begin  with  a  given  ice  pressure,  assumed  equal  to  the  overburden  pressure, 
as  in  some  cases  of  pure,  segregated  ice  near  the  top  of  the  soil.  The 
curvature  function  F(0)  then  is  assigned  a  value  in  keeping  with  a  pore 
size  of  some  "critical  radius"  characteristic  of  the  soil.  This  procedure 
yields  a  maximum  pressure  difference  and  a  value  of  IV.  via  equation 

(2) .  Given  additional  assumptions  about  the  temperature  and  pressure 
behavior  of  the  rest  of  the  system,  one  can  determine  whether  heave  occurs 
and  attempt  to  estimate  heaving  rate  or  maximum  heaving  pressures. 
Experience  suggests,  however,  that  this  approach  may  be  too  simple  and 
that  it  cannot  satisfactorily  account  for  observed  heaving  behavior. 

Kbopmans  and  Miller  (1966)  introduced  the  hypothesis  that  the  phase 
interface  configuration  at  a  given  unfrozen  water  content,  0,  could  be 
regarded  as  the  same  for  an  unsaturated  unfrozen  sample  and  a  saturated 
frozen  sample,  for  a  given  soil.  That  is,  the  same  function  F(0)  should 
apply  in  both  cases.  If  frozen  and  unfrozen  cases  are  compared,  and  IV 
is  kept  at  zero  in  each,  then  from  (2)  we  expect 

Pi/Pa  *  aw  (3) 

From  measurements  on  frozen  and  unfrozen  samples  of  the  same  granular 
soil,  Kbopmans  and  Miller  found  this  ratio  to  be  about  0.45,  which  is  in 
reasonable  agreement  with  estimates  of  the  right-hand  side  of  equation 

(3) ,  based  on  other  measurements  (Hesstvedt  1964).  Fbr  a  clay,  however, 
this  ratio  was  found  to  be  unity,  indicating  that  in  extremely  fine¬ 
grained  soils,  noncapillary  effects  dominate. 


For  unfrozen,  unsaturated  soils,  curves  representing  F(e)  are  called 
"soil  moisture  characteristic"  curves.  Values  differ  somewhat  depending 
cm  whether  a  drying  or  wetting  sequence  is  used  to  obtain  the  data 
(hysteresis  effect).  Neglecting  the  hysteresis  effect,  if  air  pressure  is 
atmospheric,  a  unique  curve  relating  moisture  content  and  pressure  in  the 
liquid  phase  is  obtained.  In  an  analogous  manner,  in  the  case  of  frozen 
saturated  soils,  one  may  construct  "freezing  characteristic"  curves  by 
relating  0  to  two  other  variables:  for  example,  either  to  and  Pw 
as  in  (2),  or  to  Pw  and  T  by  substituting  (2)  into  (1)  (assuming  that  n 
is  either  known  or  negligible).  The  phase  composition  curves  illustrated 
in  Figure  3,  relating  the  unfrozen  water  content,  Wu,  and  temperature, 
are  special  cases  of  freezing  characteristic  curves.  That  0  should  depend 
on  both  pressure  and  temperature  corresponds  to  common  experience  with 
bulk  water.  At  a  given  temperature,  pressure  melting  may  be  induced  by  an 
increase  in  applied  pressure.  Similarly,  at  a  given  pressure,  one  may 
induce  freezing  or  thawing  by  changing  the  temperature.  In  a  heaving 
soil,  both  pressure  and  temperature  are  variable,  and,  in  a  complete 
description,  both  of  these  variables  must  be  related  to  0,  simulta- 
eously. 

Misunderstanding  of  the  freezing  characteristic  curves  is  widespread 
For  example,  liquid  pressure  (gage)  in  a  freezing  soil  cannot  be  obtained 
directly  from  an  (unfrozen)  moisture  characteristic  curve.  Although 
research  indicates  that  F(o)  may  be  the  same  in  both  cases,  the  ice-water 
surface  tension,  u»iw  must  be  used  in  (2)  in  place  of  u>aw.  Moreover, 
the  quantity  wiWF(0)  yields  values  for  the  difference  between  the  ice 
and  liquid  pressures;  absolute  values  of  each  are  arbitrary.  In  some 
modeling  approaches,  it  is  assumed  that  the  ice  pressure  is  zero  (gage). 
For  example,  when  overburden  pressure  is  negligible,  it  may  be  thought 
that  IV  can  be  determined  from  (2)  with  P*  ■  0.  This  may  not  be  the 
case;  if  the  unfrozen  water  is  in  tension,  substantial  compressive  stress 
may  still  have  been  built  up  in  the  ice,  in  accordance  with  the  capillary 
pressure  jump  between  the  phases  and  the  two  pressures  may  easily  combine 
to  produce  very  low  neutral  stress.  In  principle,  an  infinite  number  of 
combinations  of  liquid  and  ice  pressure  are  possible  for  any  value  of  0. 

The  ability  to  specify  pressures  is  extremely  important  in  heave 
models  because  of  the  manner  in  which  the  flow  of  water  must  be  handled. 
The  flow  equation  is  constructed  upon  the  assumption  of  the  Darcy  law  for 
fluid  flow  and  the  conservation  of  mass.  For  saturated  unfrozen  soil,  the 
water  content,  @,  is  equal  to  the  porosity.  Hence,  it  changes  with  time 
only  if  the  soil  matrix  deforms.  Given  appropriate  boundary  and  initial 
conditions  together  with  the  parameter  values,  the  resulting  flow  equation 
can  be  solved. 

For  unsaturated  soil,  the  water  content,  0,  is  not  equal  to  the 
porosity.  Therefore,  an  appropriate  relationship  is  required  to  relate 
the  water  content  and  the  pressure  in  the  water  phase,  Pa.  One  may 
either  express  0  in  terms  of  (Pw  -  Pa) ,  where  Pa  is  atmospheric 
pressure,  or  one  may  express  IV  terms  of  0  and  Pa.  In  the  unfrozen 


unsaturated  case,  Pa  is  essentially  constant ,  and  a  unique  relation 
between  P*  and  0  results.  From  this  result  one  may  construct  a  soil 
moisture  diffuslvity  parameter.  In  several  frost  heave  models,  soil 
moisture  diffuslvity  is  used  in  the  same  manner  as  is  customary  for  frozen 
soil.  However,  the  analogy  is  not  exact.  Unlike  Pa,  the  ice  pressure 
is  not  uniform  but  varies  from  small  values  where  the  ice  content  is  very 
low  to  higher  values,  which  approach  overburden  pressures,  at  locations 
where  ice  segregation  is  occurring. 

Expressions  of  soil-water  flux-drive  relationships  that  are  more 
general  than  Darcy's  law,  have  been  developed  by  Kay  and  Groenevelt  (1974) 
and  Groenevelt  and  Kay  (1974).  The  methods  of  irreversible  thermodynamics 
are  used  to  provide  general  expressions  for  mass  fluxes  in  terms  of  both 
temperature  and  pressure  gradients,  each  multiplied  by  phenomenological 
transport  coefficients.  Additional  vigorous  thermodynamic  analyses  are 
brought  to  bear  to  express  relationships  between  the  phases,  the  nature  of 
the  transport  coefficients,  and  the  coupling  of  transport  processes.  Such 
detail  and  rigor,  however,  are  not  features  of  existing  frost  heave 
models. 

The  process  of  regelation  also  is  involved.  It  has  been  observed 
that  applied  pressures  or  temperature  gradients  may  cause  relative  motion 
between  ice  and  small  particles  completely  enveloped  in  the  ice  (Gilpin 
1979,  1980;  Romkens  and  Miller  1973).  The  processes  involved  are  not 
completely  understood,  but  it  is  clear  that  ice  must  melt  on  one  side  of 
each  enveloped  particle  and  the  resulting  liquid  be  transported  around  the 
sides  of  the  particle  to  a  location  where  it  refreezes  (Hoekstra  and 
Miller  1965;  Hoekstra  1966;  Anderson  and  Morgenstern  1973).  It  has  been 
suggested  that  this  process  operates  in  heaving  soils,  so  that  net 
movement  of  ice,  as  well  as  movement  of  the  unfrozen  water,  may  occur 
(Miller  1978).  That  regelation  occurs  gives  rise  to  a  nuntoer  of 
questions.  For  example,  when  a  pressure  gradient  is  imposed  on  a  frozen 
soil  sample  to  obtain  hydraulic  conductivity  values  for  modeling  purposes, 
to  what  extent  does  movement  of  the  ice  contribute  to  the  observed  mass 
transfer?  What  other  factors  might  affect  the  ice  movement  rate  that  were 
absent  or  unknown  during  the  hydraulic  conductivity  measurement  but  may  be 
significant  in  other  situations  for  the  parameter  values  one  applied? 
Horlguchi  and  Miller  (1980)  report  difficulties  in  their  efforts  to 
distinguish  between  hydraulic  flow  and  transport  by  regelation  in 
permeameter  experiments.  Additional  concepts  or  rules  dealing  with  this 
process  may  be  needed  to  characterize  adequately  the  dynamics  of  pore  ice 
movement  and  to  relate  it  to  the  flow  of  unfrozen  water  and  heat. 

Mathematical  expressions  for  heat  flow  needed  in  frost  heaving 
simulations  begin  with  the  principle  of  energy  conservation.  Attempts 
continue  to  develop  rigorous  procedures  for  averaging  accepted  microscopic 
statements  of  the  conservation  laws  applicable  to  a  single  phase  (e.g., 
liquid  within  a  pore)  in  such  a  way  that  valid  macroscopic  relations 
result  (Gray  and  O'Neill  1976;  Hassanlzadeh  and  Gray  1979ab,  I960). 
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Frost  heave  models  usually  begin  by  assuming  a  basic  heat  balance 
relating  the  volumetric  (sensible)  heat  capacities  of  all  phases,  the  rate 
of  change  of  sensible  heat  per  unit  volume,  the  divergence  of  the  flux  of 
sensible  heat  and  the  latent  heat  effects.  Examination  of  enthalpy 
equations  pertaining  to  simple,  single-phase  systems  reveals  that  certain 
terms  may  or  may  not  be  significant  (e.g.,  Bird,  et  al.  1960).  Fourier's 
law  of  heat  conduction  usually  is  used  to  express  the  rate  of  heat  flow. 
The  expressions  one  uses  depend  in  part  on  assumptions  about  the  mass 
transfer  processes. 

The  first  efforts  at  modeling  coupled  heat  and  moisture  transport  in 
freezing  soils  were  by  Harlan  (1973)  and  Guymon  and  Luthin  (1974).  They 
used  continuum-deterministic  approaches,  sometimes  referred  to  as 
conceptual  or  physics-based  approaches  which  have  been  preferred  by  most 
researchers.  Generally,  mathematical  models  have  included  simultaneous 
heat  and  moisture  transport  in  a  one-dimensional  column.  Although  there 
is  almost  total  agreement  that  these  two  processes  must  be  included  in  any 
model  of  frost  heave,  there  is  considerable  uncertainty  and  disagreement 
on  how  to  model  ice  segregation. 

Guymon  et  al.  (1981ab)  have  discussed  a  systems  approach  to  the 
mathematical  modeling  of  frost  heave  shown  in  Figure  4.  The  prototype 
system  S  (e.g.,  a  laboratory  soil  column  or  field  soil  column)  is 
subject  to  excitations  (or  inputs),  x,  that  are  spatially  and  temporally 
distributed.  Spatially  and  temporally  distributed  outputs,  y,  then  are 
obtained.  Inputs  or  boundary  conditions  may  be  subfreezing  temperatures, 
water-table  location,  and  surface  surcharge  (overburden).  (Xitputs  may  be 
frost  heave,  soil  pore  pressure,  temperatures,  or  ice  content.  Because  it 
usually  is  impossible  to  measure  x  exactly,  subsystem  X  indicates  a 
model  process  to  determine  an  index,  x',  of  x,  which  has  some  error.  In 
most  cases,  x  is  lumped  in  space  while  preserving  to  the  extent  possible, 
any  temporal  low-frequency  dynamic  characteristics  of  x.  Model  M  is 
assumed  to  be  a  deterministic  model  based  on  the  continuum  assumption. 
Certain  parameters  arise  in  the  model  derivation  that  purport  to 
characterize  S  (e.g.,  thermal  conductivity  or  hydraulic  conductivity). 
Subsystem  P  indicates  this  modeling  or  sampling  process,  which  -  yields 
imperfectly  known  parameters,  pi.  Therefore,  model  outputs,  y',  will  be 
imprecise  but  may  be  compared  to  imperfect  observations  of  y  for  some 
bounded  time  period  to  determine  model  uncertainty,  e(t),  where 

e(t)  -  y'(t)  -  y(t) 

(One  considers  y  as  lumped  in  order  to  make  this  computation.)  Modeling 
uncertainty  is  arbitrarily  grouped  into  four  general  areas: 

1.  Errors  ai  due  to  the  choice  of  M,  which  includes  the  choice 
of  numerical  analog; 
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Figure  4.  Schematic  diagram  of  modeling  concept. 

2.  Errors  a2  due  to  spatial  and  temporal  discretization  and 
averaging ; 

3.  Errors  a3  due  to  boundary  conditions  (i.e. ,  choice  of  X)  and 
due  to  choice  of  initial  conditions; 

4.  Errors  04  due  to  the  selection  of  (i.e.  choice  of  P). 

The  total  model  uncertainty  is  some  function  of  the  at  errors 

e(t)  *  e(ai,  a2,  a3,  a4), 

where  the  ai  errors  nay  be  interrelated  and  e  may  be  nonstationary. 

Conceptual  mathematical  models  of  nonlinear  soil  transport  processes 
may  be  made  linear  to  evaluate  their  reliability  against  closed-form 
analytical  solutions.  Basically,  this  process  is  a  standard  analysis 
technique  to  ensure  that  there  are  no  computer  code  errors,  as  most  models 
must  be  solved  on  a  computer.  The  complexity  and  nonlinear  nature  of 
models  that  deal  with  most  soil  phenomena  require  comparison  to  prototype 
situations  to  evaluate  model  precision.  Comparison  of  model  simulations 
to  prototype  data  also  permits  evaluations  of  parameter  and  boundary- 
condition  sensitivity. 

It  might  be  argued  that  precise  knowledge  of  ice  segregation 
processes  will  eliminate  uncertainty  problems;  however,  this  result  is 
unlikely.  Significant  measurement  errors  for  many  of  the  parameters  that 
are  Incorporated  in  mathematical  models  of  frost  heave  (e.g. ,  hydraulic 
conductivity  of  soil)  are  inevitable.  Consequently,  modeling  or  analysis 
of  frost  heave  must  proceed  with  a  constant  appreciation  of  the 
variability  of  key  factors.  Ideally  future  models  will  take  that 
variability  into  account. 
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Current  Models 


The  limitations  of  all  current  models  may  be  characterized  collectively  by 
simply  stating  that  at  present  no  model  enjoys  universal  or  general 
acceptance.  Frost  heave  modeling  has  been  reviewed  by  Guymon  et  al. 
(1980),  Hopke  (1980),  and  more  recently  O'Neill  (1983).  Among  the 
simplest  heave  models  are  those  of  Arakawa  (1966),  Outcalt  (1980b),  and 
Dudek  and  Holden  (1979).  Each  of  these  presupposes  simple  heat  conduction 
above  and  below  a  narrowly  circumscribed  freezing  front.  A  heat  balance 
is  constructed  at  the  front  involving  the  divergence  of  the  sensible  heat 
flux,  the  upflow  of  unfrozen  water,  and  the  release  of  latent  heat 
associated  with  any  ice  lens  or  advance  of  the  freezing  front.  If  more 
heat  is  withdrawn  from  the  freezing  front  by  conduction  than  can  be 
accommodated  by  the  upflow  of  water,  then  the  freezing  front  in  these 
models  advances.  Otherwise,  heave  occurs.  In  the  model  of  Outcalt,  the 
flow  is  driven  by  a  pressure  gradient  depending  always  on  a  single  suction 
value  selected  as  representative  of  a  partially  frozen  condition.  A 
single,  representative  hydraulic  conductivity  value  is  also  applied. 

In  IXidek  and  Holden's  approach,  a  similarly  simple  liquid  drive  is 
constructed,  in  terras  of  the  pressure  drop  across  a  small  space  near  the 
freezing  front.  This  maxinun  liquid  pressure  drop  corresponds  to  a  unique 
freezing  point  (temperature)  depression  value,  dependent  on  the  identity 
of  the  soil.  Critical  pore  radius  concepts  are  applied  to  determine 
whether  possible  heave  will  in  fact  result,  given  a  heat  deficit  in  the 
freezing  zone.  These  models  have  the  virtue  of  simplicity.  They  are 
limited,  however,  by  not  including  much  of  the  physics  specific  to  the 
processes  of  ice  segregation.  Among  other  things,  neither  ice  nor 
overburden  pressure  effects  are  accomodated. 

Harlan  (1973)  and  Guymon  and  Luthln  (1974)  develope  models  of  one- 
dimensional  heat  and  moisture  flow  in  partially  frozen  soil  based  cm 
simultaneous  heat  and  moisture  flow  equations.  Although  somewhat 
different  ancillary  equations  are  used  for  subprocesses,  the  general 
approach  is  similar.  Harlan  has  used  finite  differences,  and  Guymon  and 
Luthln  have  used  finite  elements.  Jams  and  Norm  (1980)  used  Harlan's 
(1973)  model  to  solve  more  extensively  coupled  heat  and  moisture  flow 
equations,  applied  in  greater  detail  over  a  portion  of  freezing  soil  in 
which  coupled  transfer  and  phase  change  take  place.  In  an  unsaturated 
situation,  mass  flow  is  driven  by  a  gradient  in  unfrozen  water  content. 
This  gradient  is  expressed  in  the  same  way  in  both  the  frozen  and  unfrozen 
zones,  i.e.,  a  diffusivlty  formulation  is  used.  An  experimental  curve  is 
employed  for  relating  unfrozen  moisture  content  to  temperature.  Thermal 
and  mass  transfer  coefficients  vary  throughout  the  soil  according  to 
conditions,  with  soil  moisture  diffusivlty  values  patterned  after  data 
from  unfrozen  unsaturated  soil.  Laboratory  measurements  of  water 
redistribution  in  a  freezing  soil  are  estimated  well  if  the  calculated 
liquid  flux  is  reduced  by  an  impedance  factor  ranging  between  1  and  1000. 


Transport  parameters  are  assumed  to  vary  continuously  in  space,  because  no 
Ice  segregation  occurs  and  the  soil  remains  a  continuous  medium* 

A  amber  of  other  researchers  have  followed  essentially  similar 
lines,  modifying  ancillary  relationships  and  adding  various  features 
including  segregation  and  heave  in  open  systems  (Outcalt  1976,  1977; 
Taylor  and  Luthin  1978;  Sheppard  et  al.  1978;  Berg  et  al.  1980).  Instead 
of  soil  moisture  diffusivity,  Outcalt  has  used  an  unsaturated  hydraulic 
conductivity,  with  the  same  dependency  on  pressure  in  both  frozen  and 
unfrozen  zones.  Liquid  pressure  is  obtained  from  an  expression  in  terms 
of  the  vapor  pressure  of  ice,  and  liquid  content  is  taken  as  a  function  of 
pressure  only.  Simulations  of  hypothetical  systems  are  qualitatively 
reasonable.  Taylor  and  Luthin  compare  computations  with  the  measurements 
of  Jame  and  Norum  (1980)  and  also  those  of  Dirksen  and  Miller  (1966). 
They  too  reduce  a  liquid  diffusivity  applied  in  the  frozen  soil  by  a  large 
impedance  factor.  The  calculated  mass  content  in  the  system  is 
continually  readjusted  through  an  "R  coefficient,”  to  preserve 
consistency.  Sheppard,  Kay,  and  Loch  express  phase  change  rate  as  the  sum 
of  liquid  accumulation  rate  and  the  divergence  of  liquid  flux. 
Temperature  and  pressure  are  related  by  the  Clapeyron  equation,  with  the 
ice  pressure  assumed  uniformly  zero. 

Guymon  et  al.  (in  press)  have  continued  the  development  of  a  model 
first  described  by  Berg  et  al.  (1980)  and  Guymon  et  al.  (1980).  This 
model  solves  one-dimensional  coupled  heat  and  moisture  transport 
equations.  Phase  change  effects  are  estimated  by  a  simple  heat  balance 
bookkeeping  procedure  for  a  freezing  or  thawing  zone.  Consolidation 
effects  are  modeled  on  the  basis  of  Terzaghi's  theory.  The  model  Includes 
surcharge  and  overburden  effects  by  modeling  pore  pressures  as  ice 
segregation  is  occurring.  The  model  requires  the  calibration  of  hydraulic 
conductivity  values  by  forcing  model  output  to  fit  laboratory  test  results 
at  certain  times.  The  deterministic  model  is  cascaded  with  a 
probabilistic  model  (Guymon  et  al.  1981b)  to  account  for  parameter 
uncertainty.  Verification  of  results  against  laboratory  and  field  data 
are  reported  by  Guymon  et  al.  (1981a,  and  in  press)  and  Hrcmadka  et  al. 
(1982).  This  one-dimensional  model  has  also  been  extended  to  simple  two- 
dimensional  cases  (Guymon  and  Hrcmadka  1982). 

Within  relatively  homogeneous  laboratory  soil  columns,  lenses 
sometimes  are  often  observed  spanning  much  or  all  of  the  column  width.  In 
such  cases,  these  lenses  represent  substantial  discontinuities  in  the  soil 
and  must  Intercept  most  or  all  of  the  incoming  liquid  flow.  Such  radical 
changes  In  the  scale  of  phenomena  to  be  simulated  challenge  traditional 
soil  modeling  methods.  In  one-dimensional  modeling  a  lens  blocks  flow, 
but  in  the  three-dimensional  situation  lenses  are  discontinuous.  Thus,  as 
in  the  modeling  of  atmospheric  turbulence,  for  example,  only  three- 
dimensional  models  are  completely  realistic.  Attempts  to  parameterize 
three-dimensional  processes  in  one-  or  two-dimensional  model  space  are 
always  Implicit  in  model  building,  but  the  implications  usually  are  not 
considered  explicitly.  This  fact,  together  with  previously  mentioned 
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parameterization  problems,  appears  to  set  an  upper  boundary  on  model 
validity. 

The  most  physically  elaborate  numerical  models  are  those  of  Hopke 
(1980),  Gilpin  (1980),  and  O'Neill  and  Miller  (1982,  1984).  Hopke's  model 
relates  temperatures  and  pressures  using  the  Clapeyron  equation  and 
capillary  concepts.  Use  of  the  same  freezing  characteristic  function  in 
the  unsaturated  frozen  and  unfrozen  zone,  together  with  assumptions  about 
the  ice  pressure  behavior,  leads  evidently  to  discontinuities  that 
complicate  the  computations.  Ice  pressure  is  assumed  to  be  zero  at  the 
freezing  front  and  also  within  the  frozen  soil  when  pores  are  not 
completely  filled  with  liquid  water  and  ice.  Ice  pressure  builds  as 
freezing  progresses,  and  heave  is  allowed  to  occur  wherever  ice  pressure 
equals  the  overburden  pressure.  Overburden  is  included  in  such  a  way  that 
simulated  results  agree  qualitatively  with  observation,  in  that  increased 
overburden  reduces  heave.  Unlike  previous  models,  Hopke's  can  predict 
expulsion  of  water  from  the  freezing  zone,  even  while  heave  is  occurring. 
This  phenomenon  has  actually  been  observed  particularly  in  early  time  in  a 
freezing  sequence  in  which  phase  change  is  most  rapid.  High  overburden 
pressures  were  observed  to  amplify  this  effect.  Usually  the  liquid  flow 
reverses  later  in  the  sequence  when  freezing  slows,  ultimately  flowing 
into  the  freezing  zone.  Transport  is  simulated  through  a  frozen  fringe 
that  does  not  heave.  Predictions  based  on  estimated  soil  parameters  are 
compared  with  data  from  Penner  and  Ueda  ( 1977ab ) ,  and  are  qualitatively 
good,  although  calculated  heave  is  overestimated  for  conditions  that 
entail  a  thin  fringe.  Hopke  speculates  that  other  mechanisms  not 
accounted  for  in  the  model  must  limit  heave  under  this  circumstance. 

An  attempt  to  model  geological  time  scales  has  recently  been  carried 
out  by  Outcalt  (1982).  This  model  generated  massive  ground  ice  near  the 
base  of  the  active  layer  and  clearly  illustrated  the  effects  of  water 
migration  and  heave  in  permafrost  underlying  the  active  layer. 


RESEARCH  NEEDS 


As  outlined  in  the  preceding  section,  the  widespread  availability  of 
computers  has  stimulated  mathematical  modeling  of  ice  segregation.  It  is 
apparent  that  this  effort  has  fully  exploited  and  now  has  overreached  the 


available  field  and  laboratory  data  base.  Much  more  effort  should  be 
expended  in  obtaining  good  laboratory  and  field  data  to  sustain  the 
continued  development  and  verification  of  mathematical  models.  More 
detailed  field  studies,  including  studies  of  the  stratigraphy  of  massive 
ground  ice  phenomena  and  long-term  disturbance  experiments,  can  furnish 
valuable  information. 

There  is  widespread  agreement  that  the  models  of  ice  segregation  need 
to  include  better  established  equations  of  heat  and  moisture  flow  as  well 
as  realistic  heave  and  pressure  submodels.  Beyond  this  general  agreement, 
however,  there  appear  to  be  strongly  divergent  views  on  a  number  of 
fundamental  questions  that  are  the  essence  of  the  mathematical  modeling  of 
ice  segregation.  For  example,  the  form  and  structure  required  of  models 
are  uncertain.  There  is  no  agreement  on  what  must  be  Included  or  can  be 
omitted  and  how  perceived  subprocesses  are  computationally  linked. 
Neither  is  there  agreement  on  the  mathematical  form  required  for  certain 
processes.  For  example,  some  investigators  insist  on  a  "complete" 
thermodynamic  expression  for  freezing  characteristics;  others  are 
satisfied  with  a  constant  value  with  respect  to  temperature  to  represent 
unfrozen  water  content.  In  addition,  there  is  disagreement  on  how  to 
account  for  surcharge  and  overburden  pressures  in  freezing  zones  where  ice 
segregation  is  occurring. 

Conceptual  as  well  as  procedural  difficulties  obstructing  further 
progress  are  numerous  and  formidable.  It  is  not  certain  that  a 
deterministic  formulation  adequate  for  either  the  scientific  or 
engineering  needs  can  be  developed.  Progress  has  been  achieved,  and 
exhaustive  testing  and  verification  of  the  most  promising  formulations  to 
establish  the  degree  of  applicability  to  engineering  analysis  and  design 
is  an  obvious  current  need.  This  activity  will  require  both  large-scale 
field  tests  and  extensive  laboratory  tests  coordinated  to  as  high  a  degree 
as  possible. 

Although,  as  emphasized  in  the  Introduction,  there  is  general 
understanding  and  agreement  on  most  of  the  qualitative  aspects  of  the 
processes  involved  in  ice  segregation  and  frost  heaving,  intensive  efforts 
carried  on  at  many  research  centers  throughout  the  world  concerned  with 
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constructing  numerical  simulations  of  ice  segregation  and  frost  heaving 
have  culminated  in  a  clear  identification  of  a  number  of  areas  of 
controversy  and  difficulty.  There  follows  a  sumnary  listing  of  topics  of 
particular  significance  and  concern,  some  expressed  as  questions  with 
abbreviated  conroents. 

•  How  can  the  overburden  pressure  be  related  to  the  state  of  water 
(liquid,  ice,  vapor)  and  the  soil  matrix? 

-  The  overburden  pressure  is  regarded  as  being  uniformly  transmitted 
to  the  ice  at  least  within  ice  lenses. 

-  The  overburden  pressure  is  regarded  as  being  added  to  the  liquid 
pressure  without  partition. 

-  The  overburden  pressure  is  regarded  as  being  carried  entirely  by 
matrix  without  partition  among  the  ice  or  liquid  phases. 

-  The  overburden  pressure  is  partitioned  between  ice,  liquid,  and 
matrix.  If  so,  then  there  is  a  need  for  (a)  direct  measurements 
of  partitioning  coefficients  and  (b)  scaling  (indirect 
measurements  of  partitioning  coefficients). 

•  Do  relationships  between  overburden  pressure  and  other  parameters 
(e.g.,  liquid  pressure,  ice  pressure)  under  static  conditions  also  apply 
under  dynamic  (flux)  conditions? 

•  What  are  the  most  appropriate  measurable  parameters  for  the 
characterization  of  frost  heaving? 

•  For  purposes  of  numerical  simulation,  what  are  the  criteria  for 
ice  lens  formation? 

-  Volumetric  total  water  content  exceeds  original  porosity. 

-  Volumetric  water  content  equals  or  exceeds  porosity,  and  effective 
stress  goes  to  zero. 

-  Freezing  process  reaches  state  where  it  is  dominated  by  pores  of  a 
critical  size. 

•  What  is  the  effect  of  solute  rejection  during  the  freezing  process 
on  ice  segregation? 

•  How  should  coupled  transport  processes  (liquid  water,  ice,  vapor, 
heat)  be  expressed  in  simulations  of  frozen  soils  with  and  without  ice 
lenses? 

-  Deterministic  approach:  employ  average  values  for  transport 
coefficients  and  driving  forces. 

-  Stochastic  approach:  employ  scaling,  frequency  distributions, 
autocorrelation  theory,  etc. 

•  Continuing  debate  and  differing  opinions  characterize  the  attempt 
to  classify  frozen  soils  and  their  properties  and  Internal  states. 


Geotechnical  engineers  apply  terminology  and  parameters  suitable  to 
unfrozen  soils  and  transfer  them  to  frozen  soils;  others  consider  this 
approach  inappropriate  and  use  different  procedures. 

•  Uncertainty  exists  regarding  processes  in  continuous  long-term 
frost  heave  versus  seasonal  frost  heave,  as  they  relate  to  such  things  as 
the  long-term  stability  of  cold,  buried  structures.  In  certain  cases, 
continuous  long-term  frost  heave  is  possible  and  will  seriously  threaten 
certain  structures.  Resolution  of  this  issue  requires  (a)  improvements  in 
the  theory  of  heave  processes  and  (b)  long-term  field  and  laboratory 
studies. 

•  How  can  the  location  of  lens  formation  within  the  frozen  zone  be 
accomplished,  and  the  processes  associated  with  "secondary  heaving"  be 
better  defined? 

•  What  are  the  maximum  theoretical  frost  heaving  pressures  and  the 
maxiimm  observable  frost  heaving  pressures,  and  how  are  these  pressures 
related  (a)  to  displacement  and  to  time,  (b)  to  the  concept  of  "shut-off 
pressures",  and  (c)  to  maxiimm  observable  frost  heaving  pressures  in 
laboratory  and  field  settings? 

•  To  what  extent  does  water  and  ice  migration  in  frozen  ground 
occur,  and  how  can  this  movement  be  measured  reliably? 

•  During  frost  heaving,  what  are  the  sources  of  the  energy  involved, 
and  how  are  these  energies  partitioned? 

•  How  is  regelation  Involved  in  the  buildup  and  maintenance  or 
sustenance  of  frost  heaving  pressures?  As  ice  moves  in  the  soil,  does  it 
move  as  a  rigid  block,  passing  over  soil  grains  by  regelatlon? 

•  What  determines  and  governs  the  physical  processes  of  ice 
nucleation  in  freezing  soils? 

•  Do  soil  materials  have  characteristic  frequency  distributions  and 
autocorrelative  functions  that  are  useful  in  describing  their 
characteristics  and  variability? 

•  In  short-term  laboratory  simulations  of  long-term  ice  segregation 
and  frost  heaving,  what  is  the  best  technique  for  dealing  with  or  scaling 
time? 

•  How  is  the  added  complexity  of  solute  redistribution  dealt  with? 
More  information  on  solute  rejection  and  partitioning  is  needed. 

•  Gan  we  adequately  determine  pore  water  suction  during  freezing  and 
make  determinations  of  soil  consolidation  with  models? 

Research  designed  to  answer  these  crucial  questions  and  thus  add  to 
our  current  understanding  is  needed  to  allow  further  progress  in  numerical 
simulations  and  the  mathematical  modeling  of  ice  segregation  and  frost 
heaving  and  to  enlarge  our  current  understanding  of  the  basic  governing 
principles  and  processes. 


CONCLUSIONS  AND  RECOMMENDATIONS 


As  a  result  of  this  study,  a  number  of  reconmendations  intended  to  provide 
guidance  and  assistance  to  federal  and  state  agencies,  private  firms,  and 
consortia  that  are  active  in  funding  or  conducting  research  in  the 
scientific  and  engineering  aspects  of  ice  segregation  and  frost  heaving 
can  be  advanced.  The  information  contained  in  this  report,  together  with 
these  reconmendations,  can  help  accelerate  progress  toward  a  more 
fundamental  understanding  of  the  phenomena  associated  with  ice  segregation 
and  frost  heaving  and  can  help  to  facilitate  application  of  an  improved 
understanding  to  the  planning,  design,  and  maintenance  programs  associated 
with  buildings,  facilities,  airfields,  roads,  and  utilities,  wherever 
frost  heaving  is  a  hazard. 

The  processes  of  ice  segregation  and  frost  heaving  are  complex. 
Certain  portions  of  the  problems  Involving  coupled  and  interrelated 
processes  present  unusually  formidable  difficulties;  it  is  on  these  that 
future  effort  should  concentrate. 

1.  We  reconmend  a  continued  concentrated  research  effort  to 
understand  the  basic  processes  of  ice  segregation  and  frost  heaving. 
Among  the  most  critical  unresolved  questions  and  problems  are: 

•  What  are  the  basic  relationships  that  govern  the  coupled  transport 
processes  (liquid  water,  heat,  water  vapor,  ice)  in  ice  segregation? 

•  Hew  are  overburden  pressures  related  to  ice  segregation  and  frost 
heaving? 

•  Do  the  same  basic  relationships  between  overburden  pressure  and 
other  parameters  that  characterize  frost  heaving  apply  under  both  static 
and  dynamic  conditions? 

•  What  are  the  important  basic  principles  and  relationships 
governing  long-term,  as  distinct  from  seasonal,  frost  heaving? 

•  What  are  the  maximum  attainable  frost  heaving  pressures  under 
long-term  frost  heaving  conditions  and  how  are  they  best  controlled? 

•  What  are  the  sources  of  the  energy  involved  in  frost  heaving  and 
how  are  they  partitioned  among  the  coupled  processes  involved? 

•  Solute  redistribution  occurs  during  ice  segregation.  More 
information  is  needed  on  the  basic  principles  and  proccesses. 

•  How  is  regelatlon  Involved  in  the  buildup  and  maintenance  of 
frost  heaving  pressures? 
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•  What  determines  the  temperatures  and  governs  the  phenomena  of  ice 
nucleatlon  in  freezing  soils? 

2.  We  recommend  a  continued  intensive  research  effort  to  develop 
predictive  nunerlcal  models  that  by  simulating  ice  segregation  and  frost 
heaving  can  provide  a  means  of  consolidating  all  the  processes  and 
interactions  Involved.  This  will  continue  to  aid  in  attaining  a  better 
and  more  complete  understanding  of  frost  heaving  and  for  assessing  its 
impact  in  various  settings  and  circunstances.  Among  the  most  critical 
unresolved  questions  and  problems  pertaining  to  this  aspect  are: 

•  What  are  the  essential  and  most  appropriate  parameters  required? 

•  What  are  the  essential  criteria  for  ice  segregation  and  the 
formation  of  ice  lenses  or  enlarging  ice  masses? 

•  How  should  the  coupled  transport  processes  (heat,  liquid  water, 
water  vapor,  ice)  be  most  effectively  and  appropriately  formulated? 

•  How  can  the  locations  of  enlarging  ice  lenses  be  most  effectively 
specified  and  the  processes  of  "secondary  heaving"  be  better  defined? 

•  Must  the  partition  of  energy  among  the  various  processes  involved 
be  specifically  recognized  and  dealt  with?  If  so,  how? 

•  Are  there  useful  characteristic  frequency  distributions  and  auto¬ 
correlative  functions  that  can  be  generally  agreed  upon  and  worked  into 
the  various  approaches  to  numerical  simulations  of  frost  heaving? 

•  In  both  seasonal  and  long-term  simulations,  what  is  the  best 
technique  for  dealing  with  and  scaling  time? 

Progress  in  solving  all  of  these  difficult  and  controversial 
questions  is  crucial  to  the  attainment  of  a  complete,  rational  theory  of 
ice  segregation  and  frost  heaving.  Such  a  theory,  however,  is  needed  for 
the  improved  mathematical  simulations  and  numerical  models  capable  of 
producing  reliable  assessments  of  frost  heaving  hazards  in  a  variety  of 
locations  and  circumstances  where  engineering  design  and  construction  must 
take  place. 

3.  A  common  terminology,  clearly  relating  to  and  bridging  existing 
terminologies  in  the  various  disciplines,  is  needed  to  address  the 
phenomena  of  ice  segregation  and  frost  heaving  and  to  facilitate  a  more 
efficient  and  rational  discussion  of  the  processes,  consequences,  and 
methods  of  ameliorating  their  adverse  effects.  Therefore,  we  recommend 
that  the  several  national  scientific  and  professional  societies  involved 
give  explicit  attention  to  this  need  and  encourage  efforts  to  achieve  an 
appropriate,  comnon  terminology. 

4.  Improved  instrumentation  to  measure  the  critical  phenomena 
associated  with  ice  segregation  and  frost  heaving  is  needed.  Most  of  this 
instrumentation  must  be  developed  for  unique  uses  and  circunstances.  We 
reconmend  that  this  need  be  explicitly  recognized  by  funding  bodies  ancT 
that  the  design  and  fabrication  of  needed  instrumentation,  in  addition  to 


the  procurement  of  more  conventional  and  readily  available 
instrumentation,  be  a  prominent  part  of  the  budgets  of  ice  segregation  and 
frost  heaving  research  projects. 

5.  fe  recommend  that  field  and  laboratory  studies  of  ice 
segregation  and  frost  heaving  be  better  coordinated  in  order  to  accelerate 
progress  in  the  verification  and  improvement  of  the  various  mathematical 
models  currently  available  and  to  accelerate  the  quantification  of  the 
functions  and  parameters  of  the  various  interrelated  processes.  We 
recommend  that  this  coordination  be  a  part  of  the  protocol  of  all  major 
research  projects  dealing  with  problems  of  ice  segregation  and  frost 
heaving. 

In  recent  years  professional  societies  have  organized  several 
conferences  to  facilitate  the  exchange  of  knowledge  and  data  chi  the 
problems  of  ice  segregation  and  frost  heaving.  However,  much  more  could 
be  done  at  relatively  little  cost  during  the  planning  stages  of  major 
research  projects.  The  needed  coordination  could  be  accomplished  without 
compromising  the  essential  requirements  of  obtaining  and  managing 
proprietary  data.  Indeed,  if  attention  is  given  to  this  recommendation , 
the  value  of  the  proprietary  information  obtained  could  be  vastly 
increased. 


36 


references 


Andersland,  0.  B.  and  D.  M.  Anderson  (eds.),  1978.  Geotechnical 

Engineering  for  Cold  Regions.  New  York:  McGraw-Hill. 

Anderson,  D.  M. ,  1968.  Phase  composition  of  frozen  montmor i 1 lonite-water 
mixtures  from  heat  capacity  measurements.  Proceedings  of  the  Soil 
Science  Society  of  America,  30,  670-675. 

Anderson,  D.  M. ,  1967a.  The  interface  between  ice  and  silicate  surfaces. 
Journal  of  Colloid  and  Interface  Sciences,  25,  174-191. 

Anderson,  D.  M. ,  1967b.  Ice  nucleation  and  the  substrate-ice  interface. 
Nature,  216,  563-566. 

Anderson,  D.  M.  and  P.  Hoekstra,  1965a.  Crystallization  of  clay-absorbed 
water.  Science,  149,  318-319. 

Anderson,  D.  M.  and  P.  Hoekstra,  1965b.  Migration  of  inter  lamellar  water 
during  freezing  and  thawing  of  Wyoming  bentonite.  Proceedings  of 
the  Soil  Science  Society  of  America,  29,  298-504. 

Anderson,  D.  M.  and  N.  R.  Morgen  stem,  1973.  Physics,  chemistry  and 
mechanics  of  frozen  ground.  Proceedings  of  the  Second  International 
Conference  on  Permafrost,  fakutsk,  U.S.S.R  (pp.  257-288). 
Washington,  t).CJ.  :  National  Academy  of  Sciences. 

Anderson,  D.  M.  and  A.  R.  Tice,  1972.  Predicting  unfrozen  water  contents 
in  frozen  soils  from  surface  area  measurements.  Highway  Research 
Record,  393,  12-18. 

Anderson,  D.  M. ,  P.  J.  Williams,  G.  L.  Guymon,  and  D.  L.  Kane,  1984. 
Principles  of  soil  freezing  and  frost  heaving.  Technical  Council  on 
Cold  Regions  Engineering  Monograph,  Frost  Action  and  Its  Control. 


Arakawa,  K. »  1966.  Theoretical  studies  of  ice  segregation  in  soils 


Journal  of  Glaciology.  6,  255-260. 

Banin,  A.  and  D.  M.  Anderson,  1974.  Effects  of  salt  concentration  changes 
during  freezing  on  the  unfrozen  water  content  of  porous  materials. 
Journal  of  Water  Resources  Research,  10(1),  124-128. 

Berg,  R.  L.,  G.  Guymon  and  J.  Ingersoll,  1980.  Frost  heave  in  an 
instrumented  soil  column.  Cold  Regions  Science  and  Technology, 
3,  211-221. 

Beskow,  G.,  1935.  Soil  freezing  and  frost  heaving  with  special  attention 
to  roads  and  railroads.  Swedish  Geological  Society  Yearbook 
Series  C,  26,  375. 

Bird,  R.  B.,  W.  E.  Stewart  and  E.  N.  Lightfoot,  1960.  Transport 

Phenomena.  New  York:  John  Wiley. 

Bonnerot ,  R.  and  P.  Jamet,  1979.  A  third  order  accurate  discontinuous 
finite  element  method  for  the  one-dimensional  stefan  problem. 
Journal  of  Comparative  Physics,  32,  145-167. 

Bouyoucos,  G.  J.,  1917.  Classification  and  measurement  of  the  different 
forms  of  water  in  the  soil  by  means  of  the  dilatometer  method. 
Michigan  Agricultural  Experiment  Station  Bulletin,  36,  48  pp. 

Bouyoucos,  G.  J.  and  M.  M.  WcCool,  1916.  Further  studies  on  the  freezing 
point  lowering  of  soils.  Michigan  Agricultural  Experiment  Station 
Bulletin,  31,  51  pp. 

Chamberlain,  E.  J.  and  S.  E.  Blouin,  1978.  Denslfication  by  freezing  and 
thawing  of  fine  materials  dredged  from  waterways.  Third 
International  Conference  on  Permafrost,  Vol.  1  (pp.  622-628). 
Ottawa,  Ontario:  National  Research  Council  of  Canada. 

Crank,  J. ,  1981.  How  to  deal  with  moving  boundaries  in  thermal  problems. 
In  Lewis,  Morgan,  and  Zienkiewlcz  (eds. ) ,  Numerical  Methods  in  Heat 
Transfer  (pp.  177-200).  New  York:  John  Wiley. 

Crank,  J.  and  R.  D.  Prahle,  1973.  Melting  ice  by  the  Isotherm  migration 
method.  Bulletin  of  the  Institute  of  Mathematical  Applications, 
19,  12.  - 

Derjaginn,  B.  V.  and  N.  V.  Churaev,  1978.  The  theory  of  frost  heaving. 
Journal  of  Colloid  and  Interface  Science,  67,  391-396. 


Dirkaen,  C.  mud  R.  D.  Miller,  1966.  Closed-system  freezing  of  unsaturated 
soil.  Proceedings  of  the  Soil  Science  Society  of  America,  30, 

168-173. 

Doherty,  P.  C. ,  1970.  Hot  Pipe.  USGS  Computer  Contribution  4,  Open  file 
report.  Menlo  Park,  California:  Computer  Center  Division. 

Dubikow,  G.  I.,  1982.  Paragenesis  of  sheet  ice  and  frozen  rocks  in 
Western  Siberia  (in  Russian).  In  A.  I.  FOpov  (ed.),  Sheet  Ice  of 
the  Cryollthozone  (pp.  24-42).  Yakutsk:  An  “55SS  Inst. 

Merzlolovedeniia. 

Dudek,  S.  J.  M.  and  J.  T.  Holden,  1979.  A  theoretical  model  for  frost 
heave.  In  Lewis  and  Morgan  (eds. ) ,  Numerical  Methods  in  Thermal 
Problems  (pp.  216-229).  Swansea,  UK:  Pineridge  Pressl 

Dusinberre,  G.  M. ,  1949.  Numerical  Analysis  of  Heat  Flow.  Ne»  York: 
McGraw-Hill. 

Everett,  D.  H. ,  1961.  Thermodynamics  of  frost  damage  to  porous  solids. 
Transactions  of  the  Faraday  Society,  57,  1541-1551. 

Friplat,  J.  J. ,  A.  Jelll,  C.  Poncelet,  and  J.  Andre,  1965.  Thermodynamic 
properties  of  adsorbed  water  molecules  and  electrical  conduction  in 
montmoril  Ionites  and  silicas.  Journal  of  Physical  Chemistry, 
69,  2185. 

Gilpin,  R.  R. ,  1979.  A  model  of  the  "liquid-like”  layer  between  ice  and  a 
substrate  with  application  to  wire  regelation  and  particle 
migration.  Journal  of  Colloid  and  Interface  Science,  68(2), 
235-251. 

Giipin,  R.  R. ,  1980.  A  model  for  prediction  of  ice  lensing  and  frost 
heave  in  soils.  Journal  of  Water  Resources  Research,  16, 
918-930. 


Gold,  L.  W. ,  1957.  A  possible  force  mechanism  associated  with  the 
freezing  of  water  in  porous  materials.  Highway  Research  Board 
Bulletin,  168,  65-72. 

Gray,  W.  G.  and  K.  O'Neill,  1976.  Oh  the  general  equations  for  flow  in 
porous  media  and  their  reduction  to  Darcy's  law.  Journal  of  water 
Resources  Research,  12,  148-154. 

Groenevelt,  P.  H.  and  B.  D.  Kay,  1974.  On  the  interaction  of  water  and 
heat  transport  in  frozen  and  unfrozen  soils.  II:  The  liquid  phase. 
Proceedings  of  the  Soil  Science  Society  of  America.  38,  400-404. 


•  •. 


lum  on 


nstitute  o 


logy. 


Groenevelt,  P.  H.  and  B.  D.  Kay,  1980b.  On  pressure  distribution  and 
effective  stress  in  unsaturated  soils.  Canadian  Journal  of  Soil 
Science  (submitted). 

Guymon,  G.  and  T.  V.  Hroraadka,  1982.  A  Two-Dimensional  Numerical  Model 
of  Coupled  Heat  Moisture  Transport  in  frost  Heavi 


Report.  U.S.  Army  Research  Office. 

Guymon,  G.  and  J.  N.  Luthln,  1974.  A  coupled  heat  and  moisture  transport 
model  for  arctic  soils.  Journal  of  Water  Resources  Research, 
10(5). 

Guymon,  G. ,  T.  V.  Hromadka  and  R.  L.  Berg,  1980.  A  one-dimensional  frost 
heave  model  based  upon  simulation  of  simultaneous  heat  and  water 
flux.  Obld  Regions  Science  and  Technology,  3(2&3) ,  253-262. 

Guymon,  G.,  R.  L.  Berg,  T.  C.  Johnson,  and  T.  V.  Hrcraadka,  1981a.  Results 
from  a  mathematical  model  of  frost  heave.  Transportation  Research 
Record,  809,  2-6. 

Guymon,  G. ,  M.  E.  Harr,  R.  L.  Berg,  and  T.  V.  Hrcraadka,  1981b.  A 
probabilistic-deterministic  analysis  of  one-dimensional  Ice 
segregation  in  a  freezing  soil  column.  Cbld  Regions  Science  and 
Technology,  j>,  127-140. 

Guymon,  G. ,  R.  L.  Berg,  T.  C.  Johnson,  and  T.  V.  Hromadka,  in  press. 
Mathematical  Model  of  Frost  Heave  in  Pavements.  Hanover , 

New  Hampshire:  U.S.  Army  Cold  Regions  Research  and  Engineering 
Laboratory. 

Harlan,  R.  L. ,  1973.  Analysis  of  coupled  heat-fluid  transport  in 

partially  frozen  soil.  Journal  of  later  Resources  Research,  9, 
1314-1323.  “ 


Hassanizadeh,  M.  and  V.  G.  Gray,  1979a.  General  conservation  equations 
for  multi-phase  systems:  1.  Averaging  procedure.  Advances  in  later 
Resources,  2,  131-144. 


Hassanizadeh,  M.  and  W.  G.  Gray,  1979b.  General  conservation  equations 
for  multi-phase  systems:  2.  Mass,  moron  tun,  energy,  and  entropy 
equations.  Advances  in  Water  Resources,  2,  191-204. 


Hassanizadeh,  M.  and  W.  G.  Gray,  1980.  General  conservation  equations  for 
multi-phase  systems:  3.  Constitutive  theory  for  porous  media  flow. 
Advances  in  Water  Resources,  3,  25-38. 

Hesstvedt,  E. ,  1964.  The  Interfacial  Energy,  Ice/ Water.  Publication  56 
(pp.  1-4).  Oslo,  Norway:  Norwegian  Geotechnical  Institute. 

Hoekstra,  P. ,  1965.  Conductance  of  frozen  bentonite  suspension. 

Proceedings  of  the  Soil  Science  Society  of  America,  29,  519-522. 

Hoekstra,  P. ,  1966.  Moisture  movement  in  soils  under  temperature 

gradients  with  the  cold  side  below  freezing.  Journal  of  Water 

Resources  Research,  2,  241-250. 

Hoekstra,  P. ,  1969a.  The  Physics  and  Chemistry  of  Frozen  Soils. 

Highway  Research  Board  Special  Report  103  (pp.  78-90).  Washington , 
D.C. :  National  Academy  of  Sciences. 

Hoekstra,  P. ,  1969b.  Water  movement  and  freezing  pressures.  Proceedings 
of  the  Soil  Science  Society  of  America.  33,  512-518. 

Hoekstra,  P.  and  E.  Chamber  lain,  1964.  Electro-osmosis  in  frozen  soil. 
Nature,  203,  1406. 

Hoekstra,  P. ,  and  R.  D.  Miller,  1965.  Movement  of  Water  in  a  Film 
Between  Glass  and  Ice.  ResearcE  Report  153.  Hanover , 

Hew  Hampshire:  it.S.  Xnny  Cold  Regions  Research  and  Engineering 

Laboratory. 

Hopke,  S. ,  1980.  A  model  for  frost  heave  including  overburden.  Cbld 
Regions  Science  and  Technology,  £(2&3),  111-127. 

Horiguchl,  K.  and  R.  D.  Miller,  1960.  Experimental  studies  with  frozen 
soil  in  an  "ice  sandwich"  permeameter.  Cbld  Regions  Science  and 
Technology,  3(2fc3),  177-183. 

Rromadka,  T.  7.,  G.  Guymon  and  R.  L.  Berg,  1980.  Some  approaches  to 
modeling  phase  change  in  freezing  soils.  Cold  Regions  Science  and 
Technology.  4,  137-145. 

Hromadka,  T.  V.,  G.  Guymon,  and  R.  L.  Berg,  1962.  Sensitivity  of  a  frost 
heave  model  to  numerical  method.  Cold  Regions  Science  and 
Technology,  6,  1-10. 


Jane,  Y.  V.  and  D.  I.  Norum,  1960.  Heat  and  mass  transfer  in  a  freezing 
imsaturated  porous  Bedim.  Journal  of  later  Resources  Research, 
16,  811-619.  ; 


Kaplar,  C.  V.,  19*38.  New  experiments  to  simplify  fro6t  susceptibility 
testing  of  soils.  Highway  Research  Record,  215,  48-59. 


Karpov,  E.  G.,  1961.  New  information  about  thick  layered  deposits  of 
underground  ice  in  the  northern  reaches  of  the  Yenisei  (in  Russian). 
Section  of  Glaciology  Publication  No.  41  (pp.  67-70).  aloscow: 
Academy  of  Sciences  of  the  U.S.S.R. 


Kay,  B.  D.  and  P.  H.  Groenevelt,  1974.  On  the  interaction  of  water  and 
heat  transfer  in  frozen  soils:  I.  Basic  theory,  the  vapor  phase. 
Proceedings  of  the  Socil  Science  Society  of  America,  38,  395- 


Kbnrad,  J.  M.  and  N.  R.  Morgenstern,  1981.  The  segregation  potential  of  a 
freezing  soil.  Canadian  Geotechnical  Journal,  18,  482-491. 


Koopmans,  R.  V.  R.  and  R.  D.  Miller,  1966.  Soil  freezing  and  soil  water 
characteristic  curves.  Proceedings  of  the  Soil  Science  Society  of 
America,  30,  680-685. 


Korelsha,  M.  M. ,  A.  N.  Khimenkow  and  G.  S.  Bryksina,  1981.  On  the  origin 
of  layer  deposits  of  underground  ice  in  the  northern  part  of  Western 
Siberia  (in  Russian).  Section  of  Glaciology  Wblicatlon  No.  41 
(pp.  62-66).  Moscow:  Acadeny  of  Sciences  of  the  U.S.S.R. 


loch,  J.  P.  G.  and  B.  D.  Kay,  1978.  later  redistribution  In  partially 
frozen,  saturated  silt  under  several  temperature  gradients  and  over¬ 
burden  loads.  Proceedings  of  the  Soil  Science  Society  of  America, 
42,  400-406. 


lynch,  D.  R.  and  K.  O'Neill,  1961.  Continuously  deforming  finite  elements 
for  the  solution  of  parabolic  problems,  with  and  without  phase 
change.  International  Journal  of  Numerical  Methods  for 


Mackay,  J.  R.,  1971.  The  origin  of  massive  icy  beds  in  permafrost, 
western  arctic  coast,  Canada.  Canadian  Journal  of  Earth  Sciences. 
8,  397-422. 


Mackay,  J.  R. ,  1973a.  A  frost  tube  for  the  determination  of  freezing  in 
the  active  layer  above  permafrost.  Canadian  Geotechnical  Journal, 
10,  392-396.  - 


Mackay,  J.  R.,  1973b.  Problems  In  the  origin  of  massive  Icy  beds,  western 
artic,  Cknada.  Permafrost:  The  North  American  Contribution  to  the 
Second  International  Conference  (pik  223-228).  Washington,  b.C. : 
National  Academy  of  Sciences. 

liackay,  J.  R. ,  1983a.  Downward  water  movement  Into  frozen  ground,  Western 
Arctic  Coast,  Canada.  Canadian  Journal  of  Earth  Sciences,  20, 
120-134. 

Mackay,  J.  R. ,  1983b.  Oxygen  Isotope  variations  In  permafrost, 

Tuktoyaktuk  Peninsula  Area,  Northwest  Territories.  Current 
Research.  Riper  83-1B  (pp.  67-74).  Geological  Survey  of  Canada. 

Martynov,  G.  S.,  1956.  The  calorimetric  method  of  determining  the 

quantity  of  unfrozen  water  in  frozen  soil.  In  L.A.  Melster  (ed.), 
Data  on  the  Principles  of  the  Study  of  Frozen  Zones  on  the  Earth's 
Crust.  Issue  III.  Moscow:  Academy  of  Sciences  U.S.S.R.  (National 
Research  Council  of  Canada,  Ottawa,  Technical  Translation  1088.) 

Miller,  K. ,  1981.  Moving  finite  dlements  #11  FIAM.  Journal  of  Numerical 
Analysis.  18,  1033-1057. 

Miller,  R.  D. ,  1972.  Freezing  and  heaving  of  saturated  and  unsaturated 
soils.  Highway  Research  Record,  393,  1-11. 

Miller,  R.  D. ,  1973.  Soli  freezing  in  relation  to  pore  water  pressure 
and  temperature.  Proceedings  of  2nd  International  Conference  on 
Permafrost  (pp.  344-352).  Washington,  D.C. :  National  Academy  o T 
Sciences. 

Miller,  R.  D. ,  1977.  Lens  initiation  in  secondary  heaving.  Proceedings 
of  the  International  Symposiun  cm  Frost  Action  In  SoilsT  Lulea, 
Sweden:  University  of  Lulea. 

Miller,  R.  D. ,  1978.  Frost  heaving  in  non-col loidal  soils.  Proceedings 
of  3rd  International  Conference  on  Permafrost,  Edmonton ,  Alberta , 
Canada  (pp.  707-713).  Ottawa,  Ontario:  National  Research  Council 
of  Canada. 

Miller,  R.  D. ,  1980a.  Freezing  phenomena  In  soils  (Chapter  11).  In 
Applications  of  Soil  Physics.  New  York:  Academic  Press. 

Miller,  R.  D. ,  1980b.  The  adsorbed  film  controversy.  Cold  Regions 
Science  and  Technology,  2*  83-86. 


Miller,  R.  D. ,  J.  H.  Baker,  and  J.  H.  Kblaian,  1960.  Particle  size, 
overburden  pressure,  pore  water  pressure  and  freezing  temperature  of 


ice  lenses  in  soil.  Transactions  of  the  International  Congress  op 
Soil  Science, 1( 7 ) ,  122-129. 


Mujrmann,  R.  P. ,  D.  M.  Anderson,  and  J.  V.  Peck,  1970.  Ionic  diffusion  at 
the  ice-solid  interface.  Snow  Removal  and  ice  Control  Research, 
Highway  Research  Board  Special  Report  115  (pp.  78-86).  Washington, 
D.C. :  National  Acadeiqy  of  Sciences. 

Nakano,  7.,  R.  J.  Martin  and  M.  Smith,  1972.  Ultrasonic  velocities  of  the 
dllatatlonal  and  shear  waves  in  frozen  soils.  Journal  of  later 
Resources  Research,  8,  1024-1030.  “ 

Nersesova,  Z.  A.  and  N.  A.  Tsytovich,  1966.  Unfrozen  water  in  frozen 

soils.  Permafrost;  Proceedings  of  an  International  Conference 

(pp.  230-234).  Washington^  D.C. :  National  Academy  of  Sciences. 

O'Neill,  K. ,  1983.  The  physics  of  mathematical  frost  heave  models:  A 
review.  Cold  Regions  Science  and  Technology.  6,  275-291. 

O'Neill,  K.  and  R.  0.  Miller,  1982.  Numerical  Solutions  for  a  Rigid-Ice 
Model  of  Secondary  Frost  Heave.  cRReJL  Report  82-13.  Hanover, 

New  Hampshire:  U.S.  Army  Gold  Regions  Research  and  Engineering 

laboratory. 

O'Neill,  K.  and  R.  D.  Miller,  1984.  Exploration  of  a  Rigid  Ice  Model  of 
Frost  Heave.  CRREL  Report  313.  ffenover,  New  Hampshire:  07  Amy 
Cold  Regions  Research  and  Engineering  laboratory. 

CXitcalt,  S.  I.,  1976.  A  numerical  model  of  ice  leasing  in  freezing 

soils.  Proceedings  of  the  Second  International  Conference  on  Soil- 
later  PrdDlems  in  Cold  Regions  (pp.  63-74).  Washii«ton,  D.C.: 
American  Geophysical  Union. 

Outcalt,  S.  I.,  1977.  Nunerical  modeling  of  the  ice  leasing  processes. 
International  Symposium  on  Frost  Action  in  Soils.  Vol .  2 

(pp.  75-91).  Lulea,  Sweden:  University  of  Lul  a, . 

Outcalt ,  S.  I.,  1980.  A  simple  energy  balance  model  of  ice  segregation. 
Cold  Regions  Science  and  Technology.  3(2&3),  145-152. 

(Xitcalt,  S.  I.,  1982.  Massive  near-surface  ground  ice  in  arctic  Alaska. 
Physical  Geography,  3(2),  19-29. 

Ifctterson,  D.  E.  and  M.  W.  Smith,  1980.  The  measurement  of  unfrozen  water 
content  by  time  domain  reflectoraetry:  Laboratory  test  on  frozen 
samples.  Canadian  Geotechnical  Journal,  18,  131-144. 


Fenner,  E. ,  1977.  Fundamental  aspects  of  frost  action.  International 
Symposiian  on  Frost  Action  in  Soils.  Vol.  2  (pp.  17-28).  Lulea, 
Sweden:  University  of  Luisa. 

Pehner,  E. ,  1982.  Aspects  of  ice  lens  formation .  Proceedings  of  the 

Third  International  Synposiim  on  Ground  Freezing .  Hanover ,  New 
Hampshire  (pp.  239-245).  Ranover,  Nee  Hampshire:  U.S.  Amy  Cold 
Regions  Research  and  Engineering  Laboratory. 

Fenner,  E.  and  T.  Ueda,  1977a.  Tbe  dependence  of  frost  heaving  on  load 
application.  International  Symposium  on  Frost  Action  in  Soils, 
Vol.  1  (pp.  92-101).  Lulea,  Sweden:  University  of  Lulea. 

Fenner,  E.  and  T.  Ueda ,  1977b.  Mechanical  analogy  of  a  constant  heave 

rate.  International  Symposlim  on  Frost  Action  in  Soils,  Vol.  2. 
Lulea,  Sweden:  University  of  Lulea. 

Fenner,  E.  and  T.  Ueda,  1979.  Effects  of  temperature  and  pressure  on 

frost  heaving.  Engineering  Geology,  13(1-4),  29-39. 

Fewe,  T.  L. ,  1983.  Alpine  permafrost  in  the  contiguous  United  States. 
Arctic  and  Alpine  Research,  15(2),  145-156. 

Hampton,  V.  N.  and  R.  I.  Walcott,  1974.  Gravity  profiles  across  ice-cored 
topography.  Canadian  Journal  of  Earth  Sciences,  11,  110-122. 

Romkens,  M.  J.  M.  and  R.  0.  Miller,  1973.  Migration  of  mineral  particles 
in  ice  with  a  temperature  gradient.  Journal  of  Colloid  and 
Interface  Science,  42,  103-111. 


Sheppard,  M.  I.,  B.  D.  Kay  and  J.  P.  C.  Loch,  1978.  Development  and 
testing  of  a  computer  for  heat  and  mass  flow  in  freezing  soils. 
Proceedings  of  the  Third  International  Conference  on  Permafrost, 
Vol.  1  (ro.  75-81}.  Ottawa  Ontario:  National  Research  Council  of 


Takagi,  S. ,  1960.  Sumary  of  the  adsorption  force  of  frost  heaving. 
Cold  Regions  Science  and  Technology.  3(2&3),  233-235. 

Taylor,  G.  8.  and  J.  N.  Luthin,  1978.  A  model  for  coupled  beat  and 
moisture  transfer  during  soil  freezing.  Canadian  Geotechnical 
Journal.  15,  548-555. 

Tice,  A.  R. ,  C.  M.  Burrous  and  D.  M.  Anderson,  1978.  Phase  composition 
measurements  on  soils  at  very  high  water  contents  by  the  pulsed 
nuclear  magnetic  resonance  technique.  Transportation  Research 

Record,  675,  11-14.  - - 

Tsytovich,  N.  A.,  1945.  On  the  theory  of  the  equilibrium  state  of  water 
in  frozen  soils.  Izvestiya  Akademll  Nauk  SSSR  Seriya 
Geograf icheskaya ,  9,  5-6T 

Vignes-Adler,  II. ,  1977.  On  the  origin  of  the  water  aspiration  in  a 
freezing  dispersed  medium.  Journal  of  Colloid  and  Interface 

Science,  60,  162-171. 

Vignes,  M.  K.  and  K.  Dljkema,  1974.  A  model  for  the  freezing  of  water  in 

a  dispersed  medium.  Journal  of  Colloid  and  Interface  Science. 

49,  165-172.  - - - - 

Voller,  V.  and  M.  Cross,  1981.  Accurate  solutions  of  moving  boundary 
problems  using  the  enthalpy  method.  International  Journal  of  Heat 
and  Mass  Transfer.  24,  545-556. 

Williams,  P.  J. ,  1963.  Specific  Heats  and  Unfrozen  later  Content  of 
Frozen  Soils*  Technical  Memo  76  (19.  109-126).  Ottawa,  Ontario: 
National  Research  Council  of  Canada. 


Williams,  P.  J. ,  1964a.  Experimental  determination  of  the  apparent 
specific  heat  of  frozen  soils.  Geotechnique,  14(2),  133-142. 

Williams,  P.  J. ,  1964b.  Unfrozen  water  content  of  frozen  soils  and  soil 
moisture  suction.  Geotechnique.  14(3) .  231-246. 

Williams,  P.  J.,  1967.  Properties  and  Behavior  of  Freezing  Soils. 

Report  72  (199  pp.).  Oslo,  Norway:  Norwegian  Geotechnical 
Institution.  ' 


SFTJJTTED  BIBLIOGRAPHY 


Anderson,  D.  M. ,  1977.  General  aspects  of  the  physical  state  of  water  and 
water  movement  in  frozen  soils.  International  Symposlm  on  Frost 
Action  in  Soils.  Lulea,  Sweden:  University  of  Lulea. 

Anderson,  D.  M.  and  A.  R.  Tice,  1971.  Low  temperature  phases  of 
Interfacial  water  in  clay  water  systems.  Proceedings  of  the  Soil 
Science  Society  of  America,  35,  47-54. 

Anderson,  D.  14.  and  A.  R.  Tice,  1973.  The  unfrozen  interfacial  phase  in 
frozen  soil  water  systems.  Ecological  Studies,  4,  107-124. 

Anderson,  D.  If.  and  A.  R.  Tice,  1980.  Low  temperature  phase  changes  in 
montmoril  Ionite  and  nontronlte  at  high  water  contents  and  high  salt 
contents,  (told  Regions  Science  and  Technology,  3(2&3),  139-144. 

Anderson,  D.  M. ,  P.  P.  Low  and  P.  Hoekstra,  1958.  The  density  of  water 
absorbed  by  lithium,  sodium  and  potassium-bentonite.  Proceedings  of 
the  Soil  Science  Society  of  America,  22,  99. 

Anderson,  D.  M.,  P.  P.  Lorn  and  P.  Hoekstra,  1968.  Some  thermodynamic 
relationships  for  soils  at  or  below  the  freezing  point.  Journal  of 
later  Resources  Research,  4(2),  379-394. 


47 


Anderson,  D.  M. ,  A.  R.  Tice  and  A.  Banin,  1973a.  Prediction  of  unfrozen 
water  contents  In  frozen  soils  from  liquid  limit  determinations. 
Symposium  on  Frost  Action  on  Roads,  Proceedings  Summary 


Anderson,  D.  M. ,  A.  R.  Tice  and  H.  L.  McKim,  1973b.  The  unfrozen  water 
and  the  apparent  specific  heat  capacity  of  frozen  ground.  The 
Second  International  Permafrost  Conference  ( pp .  289-294 ) . 
Washington,  D.C.:  National  Academy  of  Sciences. 

Arvidson,  W. ,  1973.  Water  flow  induced  by  soil  freezing.  M.S.  Thesis, 
Department  of  Civil  Engineering.  Alberta  Canada:  University  of 
Alberta. 

Arvidson,  W.  and  N.  R.  Morgenstern,  1974.  Water  flow  induced  by  soil 
freezing.  Proceedings  of  the  27th  Canadian  Geotechnical  Conference 
(pp.  137-143)"!  Edmonton,  Albexta,  Canada J 

Bates,  R.  E.  and  M.  A.  Bilello,  1966.  Defining  the  Cold  Regions  of  the 
Northern  Hemisphere.  Technical  Report  1757  Hanover ,  flew 
Hampshire:  U.S.  Trmy  Cold  Regions  Research  and  Engineering 

laboratory. 

Benoit,  G.  R. ,  1973.  Effect  of  freeze-thaw  cycles  on  aggregate  stability 
and  hydraulic  conductivity  of  three  soil  aggregate  sizes. 
Proceedings  of  the  Soil  Science  Society  of  America.  37(1),  3-5. 

Berg,  R.  L. ,  1974.  Energy  Balance  on  a  Paved  Surface.  Technical  Report 
226.  Hanover,  New  Hampshire:  U.S.  Army  Cold  Regions  Research  and 
Engineering  Laboratory. 

Berg,  R.  L.  and  N.  Staith,  1973.  Encountering  massive  ground  ice  during 
road  construction  in  central  Alaska.  Permafrost:  The  North 

American  Contribution  to  the  Second  Internationa 


_ .  Washington,  D.C. :  National  Academy  of  Sciences 

Berg,  R.  L. ,  G.  Guymon  and  K.  E.  Gartner,  1977.  A  mathematical  model  to 
predict  frost  heave.  International  Symposium  on  Frost  Action  in 
Soils.  Lulea,  Sweden:  University  of  Lulea. 

Berg,  R.  L. ,  G.  Guymon  and  J.  Ingersoll,  1979.  Conference  on  soil-water 
problems  in  cold  regions.  Cold  Regions  Science  and  Technology. 
2(3)  1980.  -  - 

Berg,  R.  L. ,  G.  Guymon,  T.  C.  Johnson  and  T.  V.  Hromadka,  1980. 
Mathematical  Model  to  Correlate  Frost  Heave  of  Pavements  with 


Blera&os,  M. ,  K.  Dijkema  and  D.  A.  DeVires,  1978.  Water  movement  m 
porous  media  towards  an  ice  front.  Journal  of  Hydrology.  37, 

137-148. 

Blachere  and  R.  N.  Young,  1974.  Failure  of  capillary  theory  of  frost 
damage  as  applied  to  ceramics.  Journal  of  the  American  Ceramic 
Society,  57(5). 

Bowen,  R.  M. ,  1976.  Theory  of  mixtures.  In  A.  C.  Eringen  (ed.). 
Continuum  Physics,  Vol.  3,  Part  1.  New  York:  Academic  Press. 

Bowen,  R.  M.,  1980.  Incompressible  porous  media  models  by  use  of  the 
theory  of  mixtures.  International  Journal  of  Engineering  Science, 
18,  1129-1148. 

Bresler,  E. ,  1973a.  Simultaneous  transport  of  solutes  and  water  under 
transient  unsaturated  flow  conditions.  Journal  of  Water  Resources 
Research,  9,  975-986. 

Bresler,  E. ,  1973b.  Anion  exclusion  and  coupling  effects  in  non-steady 
transport  through  unsaturated  soils.  Journal  of  the  Soil  Science 
Society  of  America,  37,  663-669. 

Bresler,  E.  and  R.  D.  Miller,  1975.  Estimation  of  pore  blockage  induced 
by  freezing  of  unsaturated  soil.  Permafrost :  The  North  American 
Contribution  to  the  Second  International  Conference,  Yakutsk 
(pp.  161-175).  Washington,  D.C7:  National  Academy  Sciences. 

Bresler,  E. ,  D.  Russo  and  R.  D.  Miller,  1978.  Rapid  estimate  of 
unsaturated  hydraulic  conductivity  function.  Journal  of  the  Soil 
Science  Society  of  America,  42,  170.  ~ 

Brownell,  D.  H.,  S.  K.  Garg  and  J.  f.  Pritchett,  1977.  Governing 
equations  for  geothermal  reservoirs.  Journal  of  Water  Resources 
Research,  13,  929-934. 

Burt,  T.  P.  and  P.  J.  Williams,  1976.  Hydraulic  conductivity  in  frozen 
soils.  Earth  Surface  Processes,  Vol.  1  (pp.  349-360).  New 
York :  John  Wiley. 

Carlson,  J.  H.  and  S.  I.  Outcalt,  1979.  The  Integrated  Simulation  of 
Soil  Heat-later  Flow,  Ice  Segregation,  and  the  Surface  Energy 
Transfer  Regime.  Manuscript. 

Gary,  J.  W. ,  1966.  Soil  moisture  transport  due  to  thermal  gradients. 
Journal  of  the  Soil  Science  Society  of  America,  30,  428-433. 


Ckry,  J  V.  and  H.  F.  Ifcyland,  1972.  Salt  and  water  movement  in 

unsaturated  frozen  soil.  Journal  of  the  Soil  Science  Society  of 
America,  36(4),  549-555.  : 

C£ry,  J.  V.  and  S.  A.  Taylor,  1962.  Thermally  driven  liquid  -and  vapor 
phase  transfer  of  water  and  energy  in  soil.  Proceedings  of  the  Soil 
Science  Society  of  America,  26,  417-420. 

Ckry,  J.  W.  and  S.  A.  Taylor,  1979.  Wkter  and  salt  movement  in 

unsaturated  frozen  silt:  Principles  and  field  observations. 
Journal  of  the  Soil  Science  Society  of  America,  43(1),  3-8. 

Qasagrande,  A.,  1932.  Discussion  of  frost  heaving.  Proceedings  of  the 
Highway  Research  Board,  11(1),  168-173. 

Casagrande,  A.,  1938.  Effects  of  Frost  In  Soils.  Permanent 

International  Association  of  Roads.  8th  Congress,  1st  and  2nd 
Sessions  (p.  10).  The  Hague. 

Casagrande ,  A.,  1947.  Classification  and  identification  of  soils. 

Proceedings  of  the  American  Society  of  Civil  Engineers,  73(6) , 
783.  ! 


Chamberlain,  E.  J. ,  1972.  The  mechanical  behavior  of  frozen  earth 

materials  under  high  pressure.  Triaxial  test  conditions. 
Geotechnique,  22(3),  469-483. 

Chamberlain,  E.  J. ,  1975.  Freeze-Thaw  Consolidation  of  Fine-grained 

Dredged  Materials.  Hanover,  New  Hampshire:  0.S.  Army  Cold  Regions' 
Research  and  Engineering  laboratory. 

Chamberlain,  E.  J.,  1980a.  Overconsolidation  effects  of  ground  freezing. 
Proceedings  of  the  Second  International  Symposium  on  Ground 
Freezing  (pp.  325-337).  Trondheim,  Norway:  The  Norwegian  Institute 
of  Technology. 


Chamberlain,  E.  J.,  1981.  Frost  Susceptibility  of  Soil,  Review  of  Index 
Test.  Monograph  81-2^  ftanover.  New  Hampshire:  tJ.S.  Army  Cold 
Regions  Research  and  Engineering  Laboratory. 


Chamberlain,  E.  J.  and  S.  E..  Blouln,  1978b.  Engineering  properties  of 
subsea  permafrost  in  the  Prudoe  Bay  region  of  Beaufort  Sea. 

- " -  — ‘  - *J — il  Conference  on  Perms  f met- 


t  ■  *. 


•  .  I '  • 


•  i  a* 


Chamber lain,  E.  J.  and  A.  J.  Gov,  1978.  Effect  of  freezing  and  thawing  on 
the  permeability  and  structure  of  soils.  proceedings  of  the 
International  Symposiun  on  Ground  Freezing  (pp.  31-44).  8ochum, 
lest  Germany:  Ruhr-University. 

Cheng,  G. ,  1978.  Heat  transfer  in  geothermal  systems.  Advanced  Heat 

Transfer,  14,  1-106. 

Cheng,  G. ,  1982.  The  forming  process  of  thick  layered  ground  ice. 

Academia  Sinica,  25(7),  777-788. 

Colbeck,  S.  C. ,  1982.  Configuration  of  ice  in  frozen  media.  Journal  of 
the  Soil  Science  Society  of  America,  133,  116-123. 

Comini ,  G. ,  S.  Delguidice,  R.  V.  Lewis  and  0.  C.  Zienkiewicz,  1974. 
Finite  element  solution  of  nonlinear  heat  conduction  problems  with 
special  reference  to  phase  change.  International  Journal  of 
Numerical  Methods  of  Engineering,  8,  613-624. 

Conference  on  Soil-later  Problems  in  Cold  Regions  (First),  May  6-7, 
1975,  Calgary,  Alberta,  Canada.  Washington,  D.C. :  American 

Geophysical  Union. 

Conference  on  Soil-later  Problems  in  Cold  Regions  (Second),  September 
1-2,  1976,  Edmonton,  Alberta,  Canada!  Washington,  D.C. :  American 
Geophysical  Union. 

Croney,  D.  and  J.  C.  Jacobs,  1967.  The  Frost  Susceptibility  of  Soils  and 
Road  Materials.  Report  I2t907  Great  Britain!  Road  Research 
Laboratory! 

DeVries,  D.  A.,  1957.  Moisture  movement  in  porous  materials  under 
temperature  gradients.  Eos  Trans.  AGU,  38(2),  222-231. 

DeVries,  D.  A.,  1958.  Simultaneous  transfer  of  heat  and  moisture  in 

porous  media.  Eos  Trans.  AGU,  39,  909-916. 

DeVries,  D.  A.,  1966.  Thermal  properties  of  soils.  In  Vanlijk  (ed.). 

Physics  of  Plant  Environment  (pp.  210-235).  Amsterdam,  North 
Holland:  VanWljk. 

Dirksen,  C. ,  1964.  Water  movement  and  frost  heaving  in  unsaturated  soil 
without  an  external  source  of  water.  Ph.D  Dissertation,  Cornell 
University. 

Drew,  D.  A.,  1971.  Averaged  equations  for  two-phase  media.  Studies  in 
Applied  Mathematics,  50,  133-166. 


51 


Everett,  D.  H.  and  J.  If.  Haynes,  1965.  Capillary  properties  of  same  model 
pore  systems  with  special  reference  to  frost  damage.  Bulletin  RILEM 
27,  31-36. 

Federal  Power  Conmission,  1977.  Heooaroendatlon  to  the  Presidents  Alaska 
Natural  Gas  Transportation  System  (Chapter  1,  p.  1).  (Reprinted 

April  1978,  by  U.S.  Department  of  Ehergy,  Federal  Ehergy  Regulatory 
Conmission.) 

Feldnan,  G.  M. ,  1967.  Moisture  migration  and  stratified  texture  in 
freezing  soils.  Journal  of  Engineering  Physics.  13,  425-429. 

Freden,  S.  (1965).  Some  aspects  on  the  physics  of  frost  heave  in  mineral 
soils.  Surface  Chemistry,  79-90. 

Freden,  S.  and  L.  Stenberg,  1980.  Frost  heave  tests  on  tills  with  an 
apparatus  for  constant  heat  flow.  Proceedings  of  the  Second 
International  Symposium  on  Ground  Freezing  (pp.  ?60-77l ) . 
Trondheim,  Norway:  The  Norwegian  Institute  of  Technology. 

Freeze,  R.  A.,  1975.  A  stochastic-conceptual  analysis  of  one-dimensional 
groundwater  flow  in  nonuniform  homogeneous  media.  Journal  of  later 
Resources  Research,  11(5),  725-741. 

Fuchs,  M.,  G.  S.  Campbell  and  R.  I.  Papendick,  1978.  An  analysis  of 
sensible  and  latent  heat  flow  in  a  partially  frozen  unsaturated 
soil.  Journal  of  the  Soil  Science  Society  of  America,  42, 
379-385. 

Fukuda,  M. ,  A.  Orhun  and  J.  N.  Luthin,  1980.  Experimental  studies  of 
coupled  heat  and  moisture  transfer  in  soils  during  freezing.  Cold 
Regions  Science  and  Technology,  3(2&3),  223-232. 

Gardner,  W.  R. ,  1958.  Some  steady-state  solutions  of  the  unsaturated 
moisture  flow  equation  with  application  to  evaporation  from  a  water 
table.  Journal  of  the  Soil  Science  Society  of  America.  85. 
228-232.  — 

Gardner,  V.  R.  and  S.  Campbell ,  1972.  Systematic  and  random  errors  in 
dual  gamma  soil  bulk  density  and  water  content  measurement. 
Proceedings  of  the  Soil  Science  Society  of  America.  36,  393-398. 


Garg,  S.  K.  and  J.  W.  Pritchett,  1977.  On  pressure-work,  viscous 
dissipation  and  the  energy  balance  relation  for  geothermal 
reservoirs.  Advances  in  Water  Resources.  1,  41-47. 


52 


GaSkln,  P.  N.  and  G.  P.  Raymond,  1973.  Pare  size  distribution  as  a  frost 
susceptibility  criterion.  Symposium  on  Frost  Action  on 
Proceedings  Sunnary  and  Discussion,  October  1-3.  Oslo,  Norway. 

Gerdel,  R.  W. ,  1969.  Characteristics  of  the  Oold  Regions.  Report  1-A. 
Hanover,  New  Hampshire:  D.S.  Army  Gold  Regions  Research  and 

Engineering  Laboratory. 

Groenevelt,  P.  H.  and  G.  H.  Bolt,  1969.  Nbn-equi librium  thermodynamics  of 
the  soil-water  system.  Journal  of  Hydrology,  7,  358-388. 

Groenevelt,  P.  H.  and  G.  H.  Bolt,  1972.  Water  retention  in  soil. 
Journal  of  the  Soil  Science  Society  of  America.  113,  238-245. 

Groenevelt,  P.  H.  and  B.  D.  Kay,  1975.  Hydrostatics  of  frozen  soil. 

Proceedings  of  the  Conference  on  Soil-Water  Problems  In  Cold 
Regions,  May  6-7,  talgary,  Alberta,  Canada  (ppl  1.93-199 ). 
Washington,  D.C. :  American  Geophysical  Union. 

Groenevelt,  P.  H.  and  B.  D.  Kay,  1978.  Applicability  of  dual  pmtna 

scanning  to  freezing  soils  and  the  problem  of  stratification. 
Journal  of  the  Soil  Science  Society  of  America,  42(6),  858-863. 

Groenevelt,  P.  H.  and  B.  D.  Kay,  1979.  Pressure  distribution  in  frozen 
soils.  Canadian  Symposium  on  Hydrology  in  Oold  Cliaates,  Vancouver 
(pp.  1-5 )~  Ottawa,  Ontario:  National  Research  Council  of  Canada^ 

Groenevelt,  P.H.  and  B.  D.  Ray,  in  press.  On  thermo-osmosis  and 
thermofiltration  in  porous  media.  Proceedings  of  10th  International 
Conference  on  Soil  Science.  Moscow,  tJ.S.S.R.  - 

Groenevelt,  P.  H.  and  J.  Parlange,  1974.  Thermodynamic  stability  of 
swelling  soils.  Soil  Science.  118.  1-5. 

Guymon,  G. ,  1975.  Soil  moisture- temperature  for  Alaskan  lowlands. 

Journal  of  the  American  Society  of  Civil  Engineers,  Division  of 
Irrigation  and  Drainage, 101 (IR3).  - 

Guymon,  G.  and  R.  L.  Berg,  1976.  Galerkln  finite  element  analog  of  frost 
heave.  Second  International  Conference  on  Soil-Water  Problems  in 
Cold  Regions.  Edmonton.  Alberta.  Canada.  Washington,  157C7T 
American  Geophysical  Union. 


Harr,  M.  E.,  1977.  Mechanics  of  Particulate  Media  (A  Probabilistic 
Approach).  New  York:  UcGraw-Hill.  - 


Hass,  V.  M. ,  1962.  Frost  action  theories  compared  with  field 

observations.  Highway  Research  Board  Bulletin,  331.  81-97.  97. 

Hill,  D.  W. ,  1977.  The  influence  of  temperature  and  load  on  moisture 
transfer  in  freezing  soil.  M.S.  Thesis,  Department  of  Civil 
Engineering,  University  of  Alberta,  Alberta,  Canada. 

Hill,  D.  W.  and  N.  R.  Morgenstern,  1977.  Influence  of  load  and  heat 
extraction  on  moisture  transfer  in  freezing  soils.  International 
Symposium  on  Fro6t  Action  in  Soils  Vol.  I,  (pp.  7 6-51)7  Lulea, 

Sweden:  University  of  Lulea. 


Hoekstra,  P. ,  1959.  Frost  Heaving  Pressures.  Research  Report  51. 
Hanover,  New  Hampshire:  U.S.  Army  Cold  Regions  Research  and 

Engineering  laboratory. 


Hoekstra,  P. ,  1964.  Thermo-electric  cooling  for  frost  effect  tests. 
Proceedings  of  the  Soil  Science  Society  of  America,  28,  716. 

Hoekstra,  P.  and  R.  D.  Miller,  1965.  Movement  of  later  in  a  Film  Between 
Glass  and  Ice.  Research  Report  153.  Hanover,  New  Hampshire:  U.S. 
Army  Cold  Regions  Research  and  Engineering  Laboratory. 


Hoekstra,  P. ,  E.  Chamberlain  and  A.  Frate,  1965.  Frost  heaving 
pressures.  Highway  Research  Board  Record,  101,  28-38. 

Hoekstra,  P. ,  P.  F.  Low  and  D.  M.  Anderson,  1968a.  Some  thermodynamic 
relationships  for  soil  at  or  below  freezing  point:  1.  Freezir*?  point 
depressions  and  heat  capacity.  Journal  of  Water  Resources 
Research,  4(2),  379-394. 

Hoekstra,  P. ,  P.  F.  Low  and  D.  M.  Anderson,  1968b.  Some  thermodynamic 
relationships  for  soil  at  or  below  the  freezing  point:  2.  Effects  of 
temperature  and  pressure  on  unfrozen  soil  water.  Journal  of  Water 
Resources  Research,  4(3),  541-544. 

Holden,  J.  T. ,  1979.  Heat  and  Moisture  Transfer  in  Unsaturated  Soils 
During  Freezing.  Technical  Report  7819,  Department  of  Tbeorectica 
Mechanics’!  Nottingham,  United  Kingdom:  University  of  Nottingham. 


Holden,  J.  T. ,  R.  H.  Jones  and  S.  J-M.  Dudek  (1980)  Heat  and  m«b  flow 
associated  with  a  freezing  front.  Proceedings  of  the  Second 
International  Symposium  on  Ground  Freezing,  (pp.  502-514). 
Trondheim,  Norway:  Norwegian  Institute  o i  Technology. 


Horlguchi,  K. ,  1975.  Relations  between  the  heave  amount  and  the  specific 
surface  area  of  powdered  materials.  Low  Temperature  Science  Series 
A,  33,  237-242.  - 


54 


Horiguchi,  K. ,  1979.  Effects  of  the  rate  of  heat  removal  on  the  rate  of 
frost  heaving.  Engineering  Geology.  13,  63-71. 

Hwang,  C.  T.,  1976.  Predictions  and  observations  on  the  behavior  of  a 
warm  gas  pipeline  on  permafrost.  Canadian  Geotechnical  Journal , 
13,  452-480. 

Hwang,  C.  T. ,  1977.  Frost  heave  design  of  a  chilled  gas  pipeline.  30th 

Canadian  Geotechnical  Conference,  Saskatoon,  Saskatoon,  Canada. 

Hwang,  C.  T.,  D.  V.  Murray  and  W.  E.  Brooker,  1972.  A  thermal  analysis 
for  structures  on  permafrost.  Canadian  Geotechnical  Journal.  9, 
36-46. 

International  Conference  on  Permafrost  (Second),  July  13-28,  1973, 

Vakutsk,  U.S.S.ft.  Washington,  D.C. :  National  Academy  o f  Sciences. 

International  Conference  cm  Permafrost  (Third),  July  10-13,  1978, 

Edmonton,  Alberta,  Vols.  1-5.  <5ttawa,  Ontario:  National  Research 
Council  of  Canada. 

International  Symposinn  on  Frost  Action  in  Soils,  Feb.  16-18,  1977, 

Vols.  1-2.  Lulea,  Sweden:  University  of  Lulea. 

International  Symposia  on  Ground  Freezing: 

(First)  March  8-10,  1978.  Bochun,  West  Germany:  Ruhr  University. 
(Second)  June  24-26,  1980,  Preprints  Vol.  1  and  Proceedings 
Vol.  2.  Trondheim,  Norway:  Norwegian  Institute  of  Technology. 
(Selected  papers  in  Engineering  Geology,  18(1-4).) 

(Third)  June  22-24,  1982.  Hanover,  New  Hampshire:  U.S.  Army  Cold 
Regions  Research  and  Engineering  Laboratory. 

Jackson,  K.  A.  and  D.  R.  Uhlanm,  1966.  Frost  heave  in  soils:  influence 
of  particles  on  solidification.  Proceedings  of  the  International 
Conference  on  Low  Temperature  Science,  from  Physics  of  Snow  and 
Ice.  ' 

Jackson,  K.  A.,  B.  Chalmers  and  D.  R.  Uhlanm,  1956.  Study  of  Ice 
Formation  in  Soils.  Technical  Report  65.  Arctic  Construction  and 
Frost  Effects  Laboratory. 

Jackson,  K.  A.,  B.  Chalmers  and  D.  R.  Uhlanm,  1958.  Freezing  of  liquids 
in  porous  media  with  special  reference  to  frost  heave  in  soils. 
Journal  of  Applied  Physics.  29(8) ,  1178-1181. 


55 


Jackson,  K.  A.,  B.  Chalmers  and  D.  R.  Dhlamn,  1964.  The  Interaction 

between  particles  and  a  solid-liquid  Interface.  Journal  of  Applied 
Physics,  35,  2986-2993. 

Jackson,  K.  A.,  B.  Chalmers  and  D.  R.  Uhlaum,  1966a.  Article  sorting  and 
stone  migration  due  to  frost  heave.  Science,  152,  545-546. 

Jackson,  K.  A.,  B.  Chalmers  and  D.  R.  Uhl  arm,  1966b.  Frost  heave  in 

soils.  Journal  of  Applied  Physics,  37,  848-852. 


Jackson,  K.  A.,  B.  Chalmers  and  D.  R.  Uhlanm,  1970. 

Theorectical  Studies  of  the  Mechanism  of  Frost  Hea 


port  199  (p.  27).  Hanover,  New  Hampshire:  U 
Research  and  Engineering  Laboratory. 


rlmental  and 


Cold  Regions 


Jacobs,  J.  C. ,  1965.  The  Road  Research  laboratory  Frost  Heave  Test. 
Great  Britain:  Department  of  Science  and  Industrial  Research,  Road 
Research  laboratory. 

Jacobs,  J.  C.  and  G.  Vest,  1966.  Investigation  into  the  effect  of 
freezing  on  a  typical  road  structure.  Great  Britain:  Road  Research 
Laboratory. 

Johnsen,  O. ,  1977.  Frost  penetration  and  ice  accumulation  in  soils. 
International  Symposium  on  Frost  Action  in  Soils,  Vol.  1 


(pp.  102-111).  Lulea,  Sweden:  university  o 

Johnson,  B.  D.  and  R.  J.  Kettle,  1980.  Frost  heaving  and  hydraulic 
conductivity.  Proceedings  of  the  Second  International  Symposium  on 
Ground  Freezing  (dp.  735-747).  Trondheim, 


rweglan  Institute 


Jones,  R.  H.  and  Berry,  1979.  Influence  of  subgrade  properties  on  frost 
heave.  Highway  Mt>llc  Works  (July),  17-2. 


Jones,  R.  H.  and  S.  J-M.  Dudek,  1979.  A  precise  cell  compared  with  other 
facilities  for  frost  heave  testing.  Transportation  Research  Board 
Record,  705,  63-72. 

Jones,  R.  H.  and  K.  G.  Hurt,  1975.  Improving  the  repeatability' of  frost 
heave  test.  Highway  and  Road  Const ruction,  43,  (1781/8),  8-13. 

Jones,  R.  H.  and  K.  G.  Hurt,  1978.  An  osmotic  method  for  determining  rock 
and  aggregate  suction  characteristics  with  application  to  frost  heave 
studies.  Journal  of  Engineering  Geology,  11,  245-252. 

Kane,  0.  L. ,  1980.  Snowmelt  Infiltration  into  Seasonally  Frozen  Soils. 
Gold  Regions  Science  and  Technology,  3(2&3) ,  153-161. 

Kane,  D.  L.  and  J.  Stein,  1983a.  Water  movement  into  seasonally  frozen 
soils.  Water  Resources  Research,  19(6),  1547-1557. 

Kane,  D.  L.  and  J.  Stein,  1983b.  Field  evidence  of  groundwater  recharge 
in  interior  Alaska.  Proceedings  of  the  International  Permafrost 

Conference,  Fairbanks,  Alaska  (pp 572-577).  “f alls  i ‘ngton ,  D.  6 . : 

National  Academy  of  Sciences . 

Kane,  D.  L.  and  J.  Stein,  1983c.  Physics  of  snowmelt  infiltration  into 
seasonally  frozen  soils.  Proceedings  of  the  American  Society  of 
Agricultural  Engineers,  Advances  in  infiltration  (ppl  178-187). 

Chicago,  Illinois:  American  Society  of  Agricultural  Engineers. 

Kaplar,  C.  W. ,  1965a.  Some  Experiments  to  Measure  the  Frost  Heaving 

Force  in  a  Silt  Sol~n  technical  Notes.  Hanover,  New  Hampshire: 

ti.S.  Army  Cold  Regions  Research  and  Engineering  laboratory. 

Kaplar,  C.  W. ,  1965b.  Stone  migration  by  freezing  of  soils.  Science, 
149(3691),  1520-1521.  ~ 

Kaplar,  C.  W. ,  1965c.  A  laboratory  Freezing  Test  to  Determine  the 
Relative  Frost  Susceptibility  of  Soils.  Technical  Notes.  Hanover, 
New  Hampshire:  tJ.S.  Army  Cold  Regions  Research  and  Ehginering 

Laboratory. 

Kaplar,  C.  W. ,  1970.  Phenomenon  and  mechanism  of  frost  heaving.  Highway 
Research  Record,  304,  1-13. 

Kaplar,  C.  W. ,  1974.  Freezing  Test  for  Evaluating  Relative  Frost 

Susceptibility  of  Various  Soils.  Technical  Report  250.  Hanover  New 
Hampshire:  t).S.  Army  Cold  Regions  Research  and  Engineering 

T  -« - _ 


•  • 


Kaplar,  C.  W.,  T.  C.  Johnson  and  R.  L.  Berg,  1975.  Roadway  Design  in 
Seasonal  Frost  Areas.  Technical  Report  259  (p.  164).  Hanover,  New 
Hampshire:  U.  §.  Army  Cbld  Regions  Research  and  Engineering 

laboratory. 

Kaplar,  C.  V.  and  J.  F.  Haley,  1952.  Cold  room  studies  of  frost 'action  in 
soils.  Highway  Research  Board,  2,  246-247. 

Kaplar,  C.  f.  and  T.  W.  Lambe,  1971.  Additives  for  modifying  the  frost 
susceptibility  of  soils.  Technical  Report  123.  Hanover  New 

Hampshire:  U.S.  Army Cold  Regions  Research  and  Engineering 

laboratory. 

Kaplar,  C.  W.  and  K.  A.  Llnell,  1968.  Lens  Formation  and  Heaving  During 
a  Freezing  of  a  Slit  Soil  Containing  Benzene^  Technical  NoteT 

Hanover,  New  Hampshire:  U.&  Army  Co  13  Regions  Research  and 

Engineering  laboratory. 

Kay,  B.  D. ,  M.  I.  Sheppard  and  J.  P.  G.  Loch,  1977.  A  preliminary 
comparison  of  simulated  and  observed  water  redistribution  in  soils 
freezing  under  laboratory  and  field  conditions.  Proceedings  of  the 
International  Conference  on  Fro6t  Action  In  SoiTs^  (ppT  29-0)7 
Lulea,  Sweden:  University  of  Lulea. 

Kerstens,  M..  1949.  Thermal  Properties  of  Solis.  Bulletin  No.  28. 

Minnesota:  Engineering  Experiment  Station. 

Kinosita,  S. ,  19 62.  Heave  force  of  frozen  soli.  Low  Temperature 

Science,  (A21). 

Kinosita,  S. ,  1975.  Soil  water  movement  and  beat  flux  In  freezing 

ground.  Proceedings  of  the  First  Conference  on  Soil-Water  Problems 
In  Cold  Regions  (pp.  33-41).  Calgary.  Alberta.  Canada.  Washington, 
D.C. :  American  Geophysical  Union. 


Kinosita,  S.,  1979.  Effects  of  Initial  soil-water  conditions  on  frost 
heaving  characteristics.  Engineering  Geology,  13(1-4). 

Kinosita,  S.  and  T.  Ohno,  1963.  Heaving  force  of  frozen  ground.  Low 

Temperature  Science.  I(A21),  117-139. 

Kinosita,  S.  and  T.  Ohno,  1966a.  Heaving  force  of  frozen  ground.  Low 

Temperature  Science,  II(A24),  285-297. 

Kinosita,  S.  andd  T.  Ohno,  1966b.  Thin  section  of  frozen  soil.  Low 
Temperature  Science,  (A24),  311-314. 


Kiooeita,  S.  and  I.  TWcesfai,  I960.  Freezing  point  depression  In  moist 
soli.  Proceedings  of  the  Second  International  Symposium  on  Ground 
Freezing  dpp.  640-45171  Trondheim,  Norway:  Tbe  Norwegian  Institute 
of  Technology. 

Kinoeita ,  S. ,  K.  Horiguchi ,  K.  Tanuma  and  T.  Ohno,  1968.  Frost  heave  In 
Kitami  (1967-1968).  Low  Temperature  Science.  (A26),  363-381. 

Kinoeita,  S. ,  M.  Fukuda,  M.  Inoue  and  K.  Takeda,  1980.  On  the  perennially 
frozen  ground  under  a  cold  storage.  Proceedings  of  the  Second 
International  Symposium  on  Ground  Freezing  (pp.  1060-1067). 
Trondheim,  Norway:  Hie  Norwegian  Institute  of  Technology. 

Konrad,  J.  M.  and  N.  R.  Morgenstern,  1980.  A  mechanistic  theory  of  ice 
lens  formation  in  fine-grained  soils.  Canadian  Geotechnical 
Journal,  17,  473-485. 

Linell,  K.  A.  and  C.  V.  Kaplar,  1959.  The  factor  of  soil  and  material 
type  in  frost  action.  Highway  Research  Board  Bulletin,  225, 
81-126. 

Linell,  K.  A.  and  C.  W.  Kaplar,  1973a.  Risk  of  uncontrolled  flow  from 
wells  through  permafrost.  Second  International  Permafrost 
Conference  (pp.  462-468).  Washington,  D.C. :  National  Academy  oF 
Sciences. 

Linell,  K.  A.  and  C.  W.  Kaplar,  1973b.  Long-term  effects  of  vegetative 
cover  on  permafrost  stability  in  an  area  of  discontinuous 
permafrost.  Second  International  Permafrost  Conference ,  (pp.  693- 
698).  Washington,  D.C. :  National  Acadeoy  of  Science. 

Linell,  K.  A.,  F.  B.  Hennion  and  G.  F.  Lobacz,  1963.  Cbrps.  of  Engineers 
pavement  design  in  areas  of  seasonal  frost.  Highway  Research 
Record.  33,  76-128. 

Loch,  J.  P.  G. ,  1975.  Secondary  Heaving;  Experiments  and  Analysis  of 
Frost  Heaving  Pressure  in  Soils.  Ph.  D.  thesis,  department  of 
Agronomy,  Cornell  University. 

Loch,  J.  P.  G. ,  1977.  Frost  heave  mechanism  and  the  role  of  the  thermal 
regime  in  heave  experiments  on  Norwegian  silty  soils.  Norwegian 
Road  Research  laboratory,  Meddelelse,  (50),  5-14. 

Loch,  J.  P.  G. ,  1978.  Thermodynamic  equilibrium  between  ice  and  water  in 
porous  media.  Soil  Science,  126,  77-80. 


Loch,  J.  P.  G. ,  1979a.  Influence  of  the  heat  extraction  rate  on  the  ice 
segregation  rate  of  soils.  Frost  1  Jord,  20,  19-30. 

Loch,  J.  P.  G. ,  1979b.  Suggestions  for  an  improved  standard  laboratory 
test  for  frost  heave  susceptibility  of  soils.  Frost  i  Jord,  20. 
33-38. 

Loch,  j.  p.  G. ,  1980.  Frost  action  in  soils.  Proceedings  of  the  Second 
International  Symposiun  on  Ground  Freezi 


im,  Norway:  The  Norwegian  Institute  of  Technology. 

Loch,  J.  P.  G.  and  R.  D.  Mille,  1975.  Test  of  the  concepts  of  secondary 
frost  heaving.  Journal  of  the  Soil  Science  Society  of  America. 
39(6),  1036-1041.  "  - 

Logsdail ,  D.  E.  and  L.  R.  Webber,  1960.  Effect  of  frost  action  structure 
of  holdion  and  clay.  Canadian  Journal  of  Soil  Science,  39(2). 
103-106.  '  -  — 


Low,  P.  F. ,  p.  Hockstra  and  D.  M.  Anderson,  1967.  Seme  The 

Relationships  for  Soils  at  or  Below  the  Freezing  Point.  5:  Effects 


on  Unfrozen  Soil  Water. 


zzz.  Hanover,  New  Hampshire:  U.S.  Amy  Cold  Regions  Research  and 
Engineering  Laboratory. 


Lunardini,  V.  J.,  1982.  Heat  Transfer  in  Cold  Climates.  New  York:  Van 
Nostrand  Reinhold. 

Luthin,  J.  and  G.  Taylor,  1976.  Numerical  results  of  coupled  heat-mass 
flow  during  freezing  and  thawing.  Proceedings  of  the  Second 
International  Conference  on  Soil-Water 


September  1-2  (pp.  155-172),  Edmonton,  Alberta 
D.C. :  American  Geophysical  Union. 

Luthin,  J.  and  G.  Taylor,  1978.  A  model  for  coupled  heat  and  moisture 
transfer  during  soil  freezing.  nanarfian  Geotechnical  Journal, 
15,  548-555.  - - 

Mackay,  J.  R. ,  1975a.  Freezing  processes  at  the  bottom  of  permafrost, 
Tuktoyaktuk  Peninsula,  District  of  Mackenzie.  Paper  75-1A  (pp.  471- 
474).  Canadian  Geological  Survey. 

Mackay,  J.  R. ,  1975b.  Seme  resistivity  surveys  of  permafrost  thickness, 
Tuktoyaktuk  Peninsula,  N.W.T.  Riper  75-1B  (pp.  177-180).  tVnadian 
Geological  Survey. 


Mackay,  J.  R. ,  1976a.  Ice  segregation  at  depth  in  permafrost. 

flaniryical  Survey  of  Canada  Paper,  76-1  A,  287-288. 

Ifcckay,  J.  R. ,  1976b.  Cryostat ic  pressures  in  nonsorted  circles. 

Canadian  Journal  of  Earth  Sciences,  13,  880-887. 

Itickay,  J.  R. ,  1977.  Permafrost  growth  and  subpermafrost  pore  water 

expulsion,  Tuktoyaktuk  Pennisula,  District  of  Mackenzie.  Geological 
Survey  of  Canada  Paper,  77-1A,  323-326. 

Mackay,  J.  R. ,  1978.  Uplift  of  objects  by  an  upfreezlng  ice  surface. 

Canadian  Geotechnical  Journal,  15(3) ,  609-613. 

Mackay,  J.  R. ,  1979.  Plngos  of  the  Tuktoyaktuk  Peninsula  area,  Northwest 
Territories.  Geographic  Physique  et  Quaternaire,  33,  3-61. 

Mackay,  J.  R. ,  1981.  Active  layer  slope  movement  in  a  continuous 

permafrost  environment,  Gary  Island,  Northwest  Territories,  Canada. 
Canadian  Journal  of  Earth  Sciences,  18,  1686-1680. 

Mackay,  J.  R. ,  1983.  Pleistocene  ice  and  sand  wedges.  Hooper  Island, 

Northwest  Territories.  Canadian  Journal  of  Earth  Sciences,  20, 
1087-1097. 

Mackay,  J.  R.  and  R.  F.  Black,  1973.  Origin,  composition ,  and  structure 

of  perennially  frozen  ground  and  ground  ice:  A  review.  In 

Permafrost:  The  North  American  Contribution  to  the  SecontT 

International  Conference  (pp.  195-192).  'Washington  D.C. :  National 
Academy  of  Sciences.  ~ 

Mackay,  J.  R.  and  D.  K.  MacKay,  1974.  Show  cover  and  ground  temperatures, 
Garry  Island,  N.V.T.  Arctic,  27,  287-296. 

Ifackay,  J.  R. ,  J.  Ostrick,  J.  Lewis  and  D.  K.  Mackay,  1979.  Frost  heave 
at  ground  temperatures  below  0  degrees  C,  Inuvlk,  Northwest 
Territories.  Geological  Survey  of  Panada  fttper,  79-1A,  403-405. 

Mackay,  J.  R.,  V.  N.  Kbnlshev  and  A.  I.  Popov,  1980.  Geologic  controls  of 
the  origin,  characteristics  and  distribution  of  ground  ice. 

Proceedings  of  the  Third  International  Conference  on  Permafrost, 
Vol.  2.  yttawa,  Ontario:  National  Research  douncil  of  CanadaT 

Mageau,  D.  W. ,  1978.  Moisture  Migration  in  Frozen  Soil*  M.S.  Thesis. 
Department  of  Civil  Engineering,  University  of  Alberta,  Edmonton, 
Alberta,  Canada. 


61 


Mageau,  D.  W.  and  N.  Ft.  Morgenstern,  1979.  Observations  on  moisture 
migration  in  frozen  soils.  Canadian  Geotechnical  Journal,  17(1), 
54-60. 


McGaw,  R.  W. ,  1972.  Frost  heaving  versus  depth  to  water  table.  Highway 
Research  Record,  393,  45-55. 

McGaw,  R.  V.,  1973.  A  macroscopic  interpretation  of  frozen  soil  texture 
as  a  function  of  the  freezing  rate.  Symposium  on  Frost  Action  on 
Roads,  Vol •  I,  (pp.  377-385).  Oslo,  fforway:  Organization  for 

Economic  Cooperative  Development. 

McGaw,  R.  V.,  S.  I.  Outcalt  and  E.  Wg,  1978.  Thermal  properties  and 
regime  of  wet  tundra  soils  at  Barrow,  Alaska.  Proceedings  of  the 
Third  International  Conference  on  Permafrost,  VolT  1,  (pp.  47-53). 
Ottawa,  Ontario:  Rational  Research  Council  of  Canada. 

McGraw,  R.  N.  and  A.  R.  Tice,  1976.  A  simple  procedure  to  calculate  the 
volume  of  water  remaining  unfrozen  in  a  freezing  soil.  Proceedings 
of  the  Second  International  Conference  on  Soil-Water  Problems  in  Cold 
Regions.  Washington,  D.fi. :  American  Geophysical  TJnion. 

Me  Roberts,  E.  C. ,  1975a.  Some  aspects  of  a  simple  secondary  creep  model 
for  deformations  in  permafrost  slopes.  Geotechnical 

Journal,  12,  98-105. 

McRoberts,  E.  C. ,  1975b.  Field  observations  of  thawing  in  soils. 

Canadian  Geotechnical  Journal,  12. 

McRoberts,  E.  C.,  1978.  Geotechnical  Engineering,  Cold  Regions 
(Qiapter  7),  Stability  in  Cold  Regions,  (PP*  336-464).  5TT 

Anders land  and  D.  M.  Anderson  (ed.).  New  York:  Mcgraw-Hill. 

McRoberts,  E.  C.  and  N.  R.  Morgenstern,  1974.  Stability  of  thawii^j 
slopes.  Canadian  Geotechnical  Journal,  11(4),  447-469. 

McRoberts,  E.  C.  and  N.  R.  Morgenstern,  1975.  Pore  water  expulsion  during 
freezing.  Qmadlan  Geotechnical  Journal.  12,  130-141. 

McRoberts,  E.  C.  and  J.  F.  Nixon,  1973.  A  study  of  some  factors  affecting 
the  thawing  of  frozen  soils.  Canadian  Geotechnical  Journal.  12. 
159-162.  — 

McRoberts,  E.  C.  and  J.  F.  Nixon,  1974.  Reticulate  ice  veins  in 
permafrost.  Northern  Canada:  Discussion.  Canadian  Geotechnical 
Journal.  12,  159-162. 


McRoberts,  B.  C.  and  J.  P.  Nixoo  ,  1975 •  Sores-}  ^Botocholcftl  observations  on 
the  role  of  surcharge  in  soil  freezing.  Proceedings  of  the 
CbnferenraonSo  11- Water  Problems  in  Gold  Regions,  Calgary,  Canada 
(pp.  42-57).  Washington,  D.d.7  American  Geophysical  Uftion. 

m 

Meyer,  G.  H. ,  1973.  Multidimensional  Stefan  problems.  SIAM  Journal  of 
Numerical  Analysis,  10,  522-538.  ~ 

Miller,  R.  D. ,  1959a.  Discussion  of  "The  Mechanism  of  Frost  Heaving  in 
Soils,"  by  E.  Renner.  Highway  Research  Bulletin.  225,  13-16. 

Miller,  R.  D. ,  1959b.  Role  of  Electric  Double  layer  in  the  Mechanian  of 
Frost  Heaving .  Research  Report  49.  Hanover,  New  Hampshire:  U.S. 
Army  Cold  Regions  Research  and  Engineering  Laboratory. 

Miller,  R.  D. ,  1970.  Ice  sandwich  functional  semi permeable  meufcrane. 
Science,  169,  584-585. 

Miller,  R.  D. ,  1973a.  The  porous  phase  barrier  and  crystallization. 
Separation  Science,  8,  521-535. 

Miller,  R.  D.,  J.  P.  Loch  and  E.  Bresler,  1966.  Phase  equilibrium  and 
soil  freezing  in  permafrost.  Proceedings  of  the  International 
Permafrost  Conference.  Washington ,  D. C.:  National  Academy  oF 

Sciences. 

Miller,  R.  D. ,  J.  P.  Loch  and  E.  Bresler,  1975.  Transport  of  water  in  a 
frozen  permeameter.  Proceedings  of  the  Soil  Science  Society  of 
America,  39,  1029-1036. 

Morgan,  K. ,  R.  W.  Lewis  and  0.  C.  Zienkiewicz,  1978.  An  improved 
algorithm  for  heat  conduction  problems  with  phase  change. 
International  Journal  of  Numerical  Methods  of  Engineering,  12, 


Morgenstern,  N.  R.  and  J.  F.  Nixon,  1971.  One  dimensional  consolidation 
of  thawing  soils.  Canadian  Geotechnical  Journal.  10,  558-565. 

Morgenstern,  N.  R.  and  D.  E.  Pufahl,  1979.  Stabilization  of  planar 
landslides  in  permafrost.  Canadian  Geotechnical  Journal.  16m. 
737-747.  - -  — 


Morgenstern,  N.  R.  and  D.  E.  Pufahl,  1960.  Remedial  measures  for  slope 
instability  in  thawing  permafrost.  Proceedings  of  the  Second 
International  Symposium  on  Ground  Freezing,  (pp.  TOftft-llfllS. 

"  itute  of  Techno1'— • 


Morgenstern,  N.  R.  and  L.  B.  Smith,  1963.  Thaw  consolidation  tests  on 
remoulded  clay.  QnnaHinn  Geotechnical  Journal.  10(1),  25-40. 


Nielsen,  D.  R. ,  J.  V.  Biggar  and  K.  T.  Erb,  1973.  Spatial  variability  of 
field-measured  soil-water  properties.  Hilgardia,  42(7). 


Nixon,  J.  F.  and  N.  R.  Morgenstern,  1973.  The  residual  stress  in  thawing 
soils.  Canadian  Geotechnical  Journal,  10(4),  571-580. 


Cbermeiler,  S.  F. ,  1973.  Frost  heave  susceptibility  research.  Symposium 
on  Fro6t  Action  on  Roads,  Report  I.  (pp.  251-266).  Oslo,  Norway: 
Organization  for  Economic  Cooperative  Development. 


Ockendon,  J.  R.  and  V.  R.  Hodgkins,  1975.  Movli 
Heat  Flow  and  Diffusion.  Oxford:  Clarendon 


Problems  in 


O'Neill,  K.,  1982.  A  fixed  mesh  finite  element  solution  for  cartesian  2-D 
phase  change.  Journal  of  Energy  Resource  Technology.  105, 
436-441.  - 


O'Neill,  K.  and  D.  R.  Lynch,  1979.  A  finite  element  solution  for  porous 
medium  freezing  using  hermite  basis  functions  and  a  continuously 
deforming  coordinate  system.  Proceedings  of  the  First  International 
Conference  on  Nunerical  Methods  in  Thermal  Problems.  Swansea: 
Pineridge  Press. 


O'Neill,  K.  and  R.  D.  Miller,  1960.  Nmerical  solutions  for  rigid  ice 
model  of  secondary  frost  heave.  Second  International  Symposium  on 
Ground  Freezing  (pp.  656-669).  Trondheim,  Norway :  Norwegian 

Institute  of  Technology. 


O'Neill,  K.  and  G.  F.  Pinder,  1961.  A  derivation  of  the  equations  for 
transport  of  liquid  and  heat  in  three  dimensions  in  a  fractured 
porous  medium.  Advances  in  later  Resources.  4,  150-164. 


Osterkamp,  T.  E.,  1975. 
ground.  Proceedli 
Cold  Regions  (pp. 
Onion. - 


Structure  and  properties  of  ice  lenses  in  frozen 
s  of  the  Conference  on  Soil-Water  Problems  in 
9-111).  Washington,  D.C. :  American  Geophysical 


Outcalt,  S.  I.,  1969.  Weather  and  diurnal  frozen  soil  structure  at 
Charlottesville,  Virginia.  Journal  of  later  Resources  Research. 
5(6),  1377-1382. 


Oitcalt,  S.  I.,  1971.  Field  observations  of  soil  temperature  and  water 
tension  feedback  effects  on  needle  ice  nights.  Archivfuer 
Meteorolgie  und  Geophvsik,  Serle  A,  20,  43-53. 


Outcalt,  S.  I.»  1979a.  The  influence  of  the  addition  of  water  vapor 
diffusion  on  the  numerical  simulation  of  the  process  of  ice 
segregation.  Frost  i  Jord,  20,  45-57. 

Outcalt,  S.  I.,  1979b.  The  effect  of  iteration  frequency  on  a  nunerical 
model  of  near  surface  ice  segregation.  Engineering  Geology,  13, 
111-121. 

Outcalt,  S.  I.,  1980a.  A  step  function  of  ice  segregation.  Proceedings 
of  the  Second  International  Symposium  on  Ground  Freezing  (pp.  5i5- 
520).  Trondheim,  Norway:  The  Norwegian  Institute  of  Technology. 

Fenner,  E. ,  1957.  Soil  moisture  tension  and  ice  segregation.  Highway 

Research  Board  Bulletin,  168,  50-54. 

Fenner,  E. ,  1958a.  Pressures  developed  in  a  porous  granular  system  as  a 
result  of  ice  segregation.  Highway  Research  Board  Report,  40, 
191-199. 

Penner,  E. ,  1958b.  Soil  moisture  suction  -  its  importance  and 

measurement.  Proceedings  of  the  American  Society  of  Testing 
Materials,  58,  1205-1^17. 

Fenner,  E. ,  1959.  The  mechanlan  of  frost  heaving  in  soils.  Highway 

Research  Board  Bulletin,  225,  1-22. 

Fenner,  E. ,  1960.  The  importance  of  freezing  rate  in  frost  action  in 
soils.  Proceedings  of  the  American  Society  of  Testing  Materials , 
60,  1151-1165. 

Penner,  E. ,  1962.  Thermal  conductivity  of  saturated  leda  clay. 
Geotechnique,  12,  168-175. 


Fenner,  E. ,  1963.  Frost  heaving  in  soils.  Proceedings  of  the  First 
International  Conference  on  Permafrost  (pp.  197-202).  Lafayette, 
Indiana:  Firdue  University. 

Fenner,  E. ,  1966a.  The  nature  of  frost  heaving  in  soils.  Proceedings  of 
the  International  Permafrost  Conference.  Washington ,  D.C. : 
National  Academy  of  Sciences. 

Fenner,  E. ,  1966b.  Pressures  developed  during  the  undirectiooal  freezing 
of  water  saturated  porous  materials.  Proceedings  of  the 
International  Conference  on  Low  Temperature  Science,  from  Riyslcs  of 
Snow  and  Ice  (pp.  146l-l4l5). 


Penner,  E. ,  1967.  Heaving  pressures  In  soils  during  undirectlonal 

freezing.  Canadian  Geotechnical  Journal.  IV(4) ,  398-408. 

Ffenner,  E. ,  1968.  Article  size  as  a  basis  for  predicting  frost  action  in 
soils.  Soils  and  Foundations,  8(4),  21-28. 

Ftenner,  E. ,  1970a.  Frost  heaving  forces  in  leda  clay.  *  Canadian 

Geotechnical  Journal,  £(8),  8-16. 

Penner,  E. ,  1970b.  Thermal  conductivity  of  frozen  soils.  Canadian 

Journal  of  Earth  Sciences,  7,  982-987. 

Penner,  E. ,  1971.  Soil  moisture  distribution  by  ice  lensing  in  freezing 
soils.  Proceedings  of  the  17th  Annual  Meeting  of  the  Canadian 
Society  of  Soil  Science  (pp.  44-62). 

Penner,  E.,  1972.  Influence  of  freezing  rate  on  frost  heaving.  Highway 
Research,  393,  56-64. 

Penner,  E. ,  1974a.  Frost  heaving  pressures  in  particulate  materials. 
Symposliin  on  Frost  Action  on  Roads,  Oslo,  Norway,  October  1-3 
(pp.  379-385). 

Penner,  E. ,  1974b.  Uplift  forces  on  foundations  in  frost  heaving  soils. 
Canadian  Geotechnical  Journal,  11.  324-338. 

Penner,  E.  and  L.  W.  Gold,  1971.  Transfer  of  heaving  forces  by  asfreezing 
to  columns  and  foundation  walls  in  frost-susceptible  roads. 
Canadian  Geotechnical  Journal,  8,  514-526. 

Penner,  E.  and  L.  E.  Goodrich,  1980.  Location  of  segregated  ice  in  frost 
susceptible  soil.  Proceedings  of  the  Second  International  Symposium 
on  Ground  Freezing  (pp.  6*26-639).  Trondheim,  Norway!  the  Norwegian 
Institute  ot  Technology. 

Penner,  E.  and  T.  Ueda,  1978.  A  froet-susceptibility  test  and  a  basis  for 
interpreting  heaving  rates.  Proceedings  of  the  Third  International 
Conference  on  Permafrost  (pp.  ’721-727) .  Ottawa,  Ontario:  National 
Research  Council  of  Canada. 

Penner,  E.  and  T.  Alton,  1978.  Effects  of  pressure  and  temperature  on 
frost  heave.  Proceedings  of  the  International  Symposium  on  Ground 
Freezing.  Bochum,  dermany:  Ruhr-Uni versi ty . 

Perfect,  E.  and  P.  J.  Williams,  1980.  Thermally  induced  water  migration 
in  frozen  soils.  Cold  Regions  Science  and  Technology.  3(2&3), 
101-100. 


66 


I 


Philip,  J.  R.t  I960.  Thermal  fields  during  regelatlon.  Oold  Regions 
Science  and  Technology,  3(2&3),  193-204. 

Radd,  F.  J.  and  D.  H.  Oertle,  1973.  Experimental  pressure  studies  of 
•  frost  heave  mechanisms  and  the  growth-fusion  behavior  of  ice  in  soils 
and  glaciers.  Second  International  Permafrost  Conference  (pp.  377- 
384).  Washington,  D.C. :  National  Academy  of  Sciences. 

Ramachandran,  N. ,  J.  P.  Gupta  and  Y.  Jaluria,  1981.  Two-dimensional 
solidification  with  natural  convection  in  the  melt  and  convective  and 
radiative  boundary  conditions.  Journal  of  Numerical  Heat  Transfer, 
4,  469-484. 

Rampton,  V.  N.  and  J.  R.  Mackay,  1982.  Massive  ice  and  icy  sediments 
throughout  the  Tuktoyaktuk  Peninsula,  Richards  Island,  and  nearby 
areas,  District  of  Mackenzie.  Geological  Survey  of  nanaHa  Paper. 
71-21,  16.  - 

Reed,  M.  A.,  D.  W.  Lovell,  A.  G.  Altschaeffl  and  L.  E.  Wood,  1979.  Frost 
heaving  rate  reedicted  from  pore  size  distribution.  Canadian 
Geotechnical  Journal,  16,  463-472. 

Reed,  M.  A.,  C.  W.  Lovell,  A.  G.  Altschaeffl  and  L.  E.  Wood,  1977.  Frost 
heaving  rate  of  silty  soils  as  function  of  pore  distribution. 
Highway  Research  Project  Report,  77-15,  116. 

Rolph,  M.  D. ,  III  and  K.  J.  Bathe,  1982.  An  efficient  algorithm  for 
analysis  of  nonlinear  heat  transfer  with  phase  change. 
International  Journal  of  Nianerlcal  Methods  of  Engineering,  18, 
TI5P134.  ^ 

Saitoh,  T. ,  1978.  Mmerlcal  method  for  multi-dimensional  freezing 

problems  in  arbitrary  domains.  Journal  of  Numerical  Heat  Transfer, 
100,  294-299.  - 

Shamsundar,  N.  and  E.  M.  Sparrow,  1975.  Analysis  of  multi-dimensional 
phase  change  via  enthalpy  method.  Journal  of  Numerical  Heat 
Transfer,  97,  333-340. 

Sheppard,  M.  I.,  1977.  A  computer  simulation  of  freezing  soil: 
development  and  validation  using  experimental  data.  Fh.D.  Thesis 
(pp.  171),  University  of  Guelph. 

Skaven-Haug,  1956.  Frost  accumulating  foundation  mat  for  reducing  depth 
of  excavation.  Byggnadskonst ,  4(4),  91-94. 


u 


£ 

W 


V- 


Skaven-Haug ,  1959.  Protection  against  frost  heaving  on  the  Norwegian 
railroads.  Geotechnique.  IX(3),  87-106. 

Sparrow,  E.  M.,  S.  V.  Patankar  and  S.  Ramadhyani,  1977.  Analysis  of 
melting  in  the  presence  of  natural  convection  in  the  melt  region. 
Journal  of  Mechanical  Heat  Transfer.  99,  520-526. 

Stein,  J.  and  D.  L.  Kane,  1983.  In  press.  Monitoring  the  unfrozen  water 
content  of  soil  and  snow  using  time  domain  ref lectometry .  Journal 
of  Water  Resources  Research,  6(2),  1573-1584.  — — — ~ 

Stenberg,  L. ,  1980.  Frost  heave  studies  by  natural  freezing. 
Proceedings  of  the  Second  International  Symposium  on  Ground 
Freezing  (ppl  784-79471  trondnsdm,  Norway:  Norwegian  Institute  of 
technology. 

Sutherland,  H.  B.  and  P.  N.  Gaskin,  1973a.  Pore  water  and  heaving 
pressures  developed  in  partially  frozen  soils.  Second  International 
Permafrost  Conference  (pp.  409-419).  Washington,  D.C. :  National 
Academy  of  Sciences. 

Sutherland,  H.  B.  and  P.  N.  Gaskin,  1973b.  A  comparison  of  the  TRRL  and 
CRREL  tests  for  frost  susceptibility  of  soils.  Canadian 
Geotechnical  Journal,  10(3),  553-555. 

Taber,  S. ,  1930a.  Freezing  and  thawing  of  soils  as  factors  in  the 
destruction  of  road  pavements.  Publication  Roads.  VIIf6). 
113-132.  - 

Thber,  S. ,  1943.  Prennially  frozen  ground  in  Alaska.  Geological  Society 
of  America  Bulletin,  54,  1435-1548.  -  - 

Takagl,  S. ,  1963.  Theory  of  freezing-point  depression  with  special 
reference  to  soil  water.  Proceedings  of  the  Third  International 
Conference  on  Permafrost .  (pp.  216-224).  Ottawa,  Ontario:  National 
Research  Council  of  Canada. 


Takagl,  S. ,  1965.  Principles  of  Frost 
(p.  24).  Hanover ,  New  H 
and  Engineering  Laboratory. 


-luU'ciiT 


re: 


Heaving.  Report  140,  AD626174 
t).S.  Aiw  Cold  Regions  Research 


Takagl,  S. ,  1978.  Segregation  Freezing  as  the  Cause  of  Suction  force  for 
Ice  Ley  Formation.  Report  78-6.  Hanover,  New  Hampshire:  07 
Aray  Cold  Regions  Research  and  Engineering  laboratory. 

Takagl,  S. ,  1980.  The  Adsorption  Force  Theory  of  Frost  Beavl»g 

2  (PP*  57-81).  Hanover,  New  Hampshire:  Army  Cold  Regions  Research 
and  Engineering  laboratory. 


68 


t~T 


n  •  » 


il  Freezln 


of  Frozen 


ghway  Research  Board  Special 
tional  Academy  of  Sciences. 


Thkashi,  T.,  1978.  Joining  of  two  tunnel  shields  under  artificial 
underground  freezing.  Proceedings  of  the  International  Symposiun  on 
Ground  Freezing.  Rochum,  Germany:  Ruhr-tJniversity. 


Takashi,  T. ,  M.  Masuda  and  H.  Yamamoto,  1974.  Experimental  study  on  the 
influence  of  freezing  speed  upon  frost  heave  ratio  of  soil  under 
constant  stress.  Journal  of  the  Japanese  Society  of  Snow  and  Ice. 
36,  49-67.  - 


Takashi ,  T. ,  T.  Ohrai,  H.  Yamaoto,  1977.  Freezing  expansion  of  soils  in 
frozen  ground  working  method.  Soil  and  Foundation,  25. 


Takashi,  T. ,  H.  Yamamoto,  T.  Ohrai  and  M.  Masuda,  1978.  Effect  of 

penetration  rate  of  freezing  and  confining  stress  on  the  frost  heave 
ratio  of  soil.  Proceedings  of  the  Third  International  Conference  on 
Permafrost,  Vol .  F)  (pp.  736-742).  Ottawa,  Ontario:  National 

Research  Council  of  Canada. 


Takashi,  T. ,  T.  Ohrai,  H.  Yamaoto  and  J.  Okamoto,  1979.  Upper  limit  of 
heaving  pressure  obtained  by  observing  pore  water  pressure  under 
partial  soil  freezing.  SEPPTO,  Journal  of  the  Japanese  Society  of 
Snow  and  Ice,  41(4),  47-577 


Takashi,  T. ,  T.  Ohrai,  H.  Yamaoto  and  J.  Okamoto,  1980.  Upper  limit  of 
heaving  pressure  derived  by  pore  water  pressure  measurements  of 
partially  frozen  soils.  Proceedings  of  the  Second  International 
l\xn  on  Ground  Freezl 


rweglan  Institute  of  Technology. 


Thompson,  M.  R. ,  1969.  Mixture  Design  for  Lime-Treated  Soils.  Civil 
Engineering  Studies,  Highway  Engineering  Series  22. 


Thompson,  M.  R. ,  1973.  IXirability  and  frost  resistance  of  lime  and  cement 
stabilized  soils.  Symposium  on  Frost  Action  on  Roads.  Oslo, 
Norway:  Organization  ror  Economic  Cooperative  {Development. 


Thompson,  M.  R.  and  B.  J.  Dempsey,  1970.  Quantitative  Charae tarty** inn 
of  Cyclic  Freezing  and  Thawing  in  Stabilized  Pavement  Ihterials. 
Research  Record  No.  304  (pp.  33-44!).  Washington,  D.c. :  Highway 
Research  Board. 

Tice,  A.  R. ,  D.  If.  Anderson  and  A.  Banin,  1976.  The  Prediction  of 
Unfrozen  Water  Contents  in  Frozen  Soils  from  Liquid  Limit 
Deteralnat ions .  Report  76-8.  Hanover,  New  Hampshire:  U.S.  Army 
Cold  Regions  Research  and  Engineering  laboratory. 

Tice,  A.  R. ,  C.  M.  Burro  us  and  D.  M.  Anderson,  1978.  Determination  of 
unfrozen  water  in  frozen  soil  by  pulsed  nuclear  magnetic  resonance. 
Proceedings  of  the  Third  International  Conference  on  Permafrost, 
(pp.  149-155).  Ottawa,  Ontario:  National  Research  Council  of 

Canada. 

Warrick,  A.  W.  and  D.  R.  Nielsen,  1980.  Spatial  Variability  of  Soil 
Physical  Properties  in  the  Field,  in  Applications  of  Soil  Physics, 
D.  dillel  (ed.).  New  York:  Academic  Press. 

Williams,  P.  J. ,  1963a.  Specific  heats  and  unfrozen  water  content  of 
frozen  soils.  Proceedings  of  the  First  Canadian  Conference  on 
Permafrost,  Memorandum  No.  706,  (pp.  109-126).  Ottawa,  Ontario: 
National  Research  Council  of  Canada. 

Williams  P.  J. ,  1963b.  Suction  and  its  effects  in  unfrozen  water  of 
frozen  soils.  Proceedings  of  the  International  Conference. 
Washington,  D.C. :  National  Academy  of  Sciences.  1 

Williams,  P.  J. ,  1964.  Unfrozen  water  content  of  frozen  soils  and  soil 
moisture  suction.  Geotechnique .  14(3),  231-246. 

Williams,  P.  J.,  1966.  R>re  pressures  at  a  penetrating  frost  line  and 
their  prediction.  Geotechnique,  16(3),  187-208. 

Williams,  P.  J. ,  1968a.  Thermoelectric  cooling  for  precise  temperature 
control  of  frozen  and  unfrozen  soils.  Canadian  Geotechnical 
Journal.  5(4),  264-266.  - 

Williams,  P.  J.,  1968b.  Ice  distribution  on  permafrost  profiles. 
Canadian  Journal  of  Earth  Sciences.  5(12),  1381-1386. 

Williams,  P.  J. ,  1972.  Use  of  the  Ice-water  Surface  Tension  Concep 
Engineering  Practice!  High  Research  Board  Research  Report, 

(pp.  19-29).  Washington,  D.C.:  National  Academy  of  Sciences. 


it  in 


Williams,  P.  J«»  1976.  Volume  change  in  frozen  soils.  Laurits  Bjerrum 
Memorial  Volume  (pp.  233-246).  Oslo,  Norway:  Norwegian 


Geotechnical  institute. 

Williams,  P.  J.,  1977.  Thermodynamic  conditions  for  ice  accumulation  in 
freezing  soils.  International  Symposium  on  Frost  Action  in  Soils. 
Lulea,  Sweden:  University  of  Lulea. 

Williams,  P.  J. ,  1979.  Pipelines  and  Permafrost  (pp.  103).  New  York: 
Longman. 

Williams,  P.  J. ,  1980.  Design  considerations  for  large-diameter  pipelines 
in  cold  regions.  Proceedings  of  the  Second  International  Symposium 
on  Ground  Freezing,  (pp.  1068-1075).  Trondheim  Norway:  Norwegian 
Institute  of  Technology. 

Williams,  P.  J.  and  T.  Burt,  1974.  Measurement  of  hydraulic  conductivity 
of  frozen  soils.  Canadian  Geotechnical  Journal,  11,  647-650. 

Wilson,  D.  G. ,  A.  D.  Solomon  and  P.  T.  Boggs,  1978.  Moving  Boundary 
Problems.  New  York:  Academic  Press. 

Wissa,  A.  and  R.  T.  fthrtin.  1968.  Behavior  of  Soils  Under  Flexible 
Pavements — Development  of  Rapid  Frost  Susceptibility  festsT 
Research  Report  BBS-77.  Cambridge,  Mass.:  Massachusetts  Institute 
of  Technology.  (Soils,  224,  1-113. 

Yong,  R.  N. ,  1967.  Ota  the  relationship  between  partial  soil  freezing  and 
surface  process.  International  Conference  Low  Temperature  Science, 
1.(2),  1375-1385. 

Yong,  R.  N. ,  1970.  Modification  of  Frost  Heaving  with  Salt — Part  II. 
Report  No.  FH-SSI-3T 

Yong,  R.  N. ,  1971.  Heave  and  heaving  pressures  in  frozen  soils. 

Canadian  Geotechnical  Journal,  8,  272-282. 

Yong,  R.  N. ,  1973.  Some  aspects  of  surflcial  salt  treatment  of 

attenuation  of  frost  heaving.  Second  International  Permafrost 
Conference,  (pp.  426-434).  Washington,  D.C. :  National  Academy  of 
Sciences. 

Yong,  R.  N. ,  1974.  Salt  migration  and  frost  heaving  of  salt  treated  soils 
la  view  of  freezing  and  thawing.  Symposium  on  Frost  Action  on 
Roads.  Oslo,  Norway:  Organization  for  Cooperative  Economic 

Development. 


ft*  * 


te  I 


Yong,  R.  N.,  1975.  Water  and  salt  redistribution  In  freezing  soils. 
Proceedings  of  the  Conference  on  Soil-Water  Problems  In  Cold 
Regions  (pp.  58-65T:  Washington.  D.C.:  American  Geophysical  tJnion. 


Yong,  R.  N,  1976.  Field  experimentation  with  chemical  alleviation  of 
frost  damage.  Proceedings  of  the  Second  Conference  on  Soil-Water 
Problems  in  Cold  Regions  (pp.  173-185).  Washington,  D.C.:  American 
Geophysical  Union. 

Yong,  R.  N.  and  C.  H.  Cheung ,  1979.  Prediction  of  salt  influence  on 
unfrozen  water  content  in  frozen  soils.  Journal  of  Ehglneerlne 
Geology.  13,  137-155.  - * 

Yong,  R.  N.  and  N.  Serag-Eldin,  1980.  Salt-treatment  effects  on  frost 
heave  performance.  Proceedings  of  the  Second  International 
Symposium  on  Ground  Freezing  (pp.  680-691).  Trondheim '  Norway: 
Norwegian  Institute  of  Technology. 

Zhestkova,  T.  N. ,  1978.  Results  of  experimental  studies  of  the  freezing 
process  in  very  fine-grained  soils.  Proceedings  of  the  Third 

International  Conference  on  Permafrost.  Vol.  1,  (pp.  155-162). 
Ottawa,  Ontario: National  Research  Council  of  Canada. 

Zhestkova,  T.  N.  and  V.  G.  Guzhov,  1976.  Sane  experimental  data  on  the 
dynamics  of  cryogenic  textures  in  freezing  ground.  Cryogenic 
Studies,  Vol.  15.  Moscow,  D.S.S.R. :  Moscow  University  Press. 

Zoller,  H.,  1973.  Frost  heave  and  the  rapid  frost  heave  test. 
Publication  Roads,  37,  211-220. 


